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1. CURRENT TECHNICAL OBJECTIVES 


1. Optimal Utilization of Laser and VLBI Observations for Reference 
Frames for Geodynamics (Grant NSG 5265) 

2. Utilization of Range Difference Observations In Geodynamics 
(Contract NAS 5-25888) 

3. Development of Models for Ice Sheet and Crustal Deformations 
(Grant NSG 5265) 


2. ACTIVITIES 

2,1 Effects of Adopting New Precession, Nutation and Equinox Corrections 
on the Terrestrial Reference Frame 


A paper on this topic was presented at the XVII General Assembly of 
the International Astronomical Union, Patras, Greece, August 17-26, 1982, 
and appears in its entirety below. It will also appear in Bulletin 
Geodesique . 




PREFACE 

These projects are under the supervision of Professor Ivan I. 
Mueller, Department of Geodetic Science and Surveying, The Ohio State 
University. The Science Advisor of RF 711055 is Dr. David E. Smith, 
Code 9;1 , Geodynamics Branch, and the Technical Officer is Mr. Jean 
Welker, Code 903, Technology Applications Center. The Technical 
Off kar for RF 712407 is Mr. C. Stephanides, Code 942. The latter 
three are at NASA/GSFC, Greenbelt, Maryland 20771. 
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EFFECTS OF ADOPT I MG NEW PRECESSION, NUTATION AND EQUINOX CORRECTIONS 
ON THE TERRESTRIAL REFERENCE FRAME 1 

Sheng-Yuan Zhu* and Ivan I. Mueller 
Department of Geodetic Science and Surveying 
The Ohio State University, Columbus, Ohio 43210 USA 


ABSTRACT. First, the paper is devoted to the effects of adopting new 
definitive precession and equinox corrections on the terrestrial reference 
frame: The effect on polar motion is a diurnal periodic term with an amplitude 

increasing linearly in time; on UT1 it is a linear term. Second, general 
principles are given the use of which can determine the effects of small rota- 
tions (such as precession, nutation or equinox corrections) of the frame of a 
Conventional Inertial Reference System (CIS) on the frame of the Conventional 
Terrestrial Reference System (CTS). Next, seven CTS options are presented, 
one of which is necessary to accommodate such rotations (corrections). The 
last of these options requiring no changes in the origin of terrestrial 
longitudes and in UT1 is advocated; this option would be maintained by even- 
tually referencing the Greenwich Mean Sidereal Time to a fixed point on the 
equator, instead of to the mean equinox of date, the current practice. 
Accomodating possible future changes in the astronomical nutation is discussed 
in the last section. The Appendix deals with the effects of differences which 
may exist between the various CTS's and CIS's (inherent in the various observa- 
tional techniques) on earth rotation parameters (ERP) and how these differences 
can be determined. It is shown that the CTS differences can be determined 
from observations made at the same site, while the CIS differences by comparing 
the ERP's determined by the different techniques during the same time period. 

INTRODUCTION 

New general precession and equinox corrections are being introduced in the 
1984 star catalogues and ephemerides. These corrections in turn will affect 
the earth rotation parameters (ERP), i.e., polar motion coordinates and UT1, and 
thus may change the frame of the Conventional Terrestrial Reference System 

*0n leave from Shanghai Observatory, China. 

Presented at XVIII General Assembly of the International Astronomical Union, 
Patras, Greece, August 17-26, 1982. 
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(CTS) [Mueller 1981]. Williams and Melbourne [1981] have already given a 
detailed discussion of these effects on UT1 and on the origin of terrestrial 
longitudes. In fact, it was this work which gave us motivation to expand the 
discussion to include the effects on all ERP's and offer additional options 
on how the necessary changes in the CTS could be accommodated. The approach 
is strictly geometric, i.e., we try to answer the question how definitive 
corrections to precession, nutation and the equinox affect the ERP's, and thus 
the CTS. Williams and Melbourne [1981] emphasize the point of how UT1 and the 
origin of longitudes will be affected in the future by the uncertainties in the' 
newly adopted corrections or how these corrections can be improved in the future 
from ensembles of Very Long Baseline Interferometer (VLB I ) or Lunar Laser Ranging 
(LLR) observations, with the desire that no or minimum additional changes result 
in the CTS. They assume that VLBI sources are observed randomly over the sky, 
while LLR observations are equally distributed only along the ecliptic, and 
therefore the resulting equations defining the changes of the origin of terres- 
trial longitudes and UT1 are technique dependent, whereas ours are not. (Putting 
it differently, they imply that if the analysis of future VLBI or LLR ensemble 
observations indicate necessary changes in UT1 and in the origin of terrestrial 
longitudes, such changes are due to the still existing imperfections in the 
newly adopted corrections to precession, equinox, etc., and when determined 
they will be biased with respect to each other because of the different sensi- 
tivities of the two ensembles of observations.) This difference in the results 
should not confuse the reader who recognizes the different purposes for which 
these papers were written. 

1. EFFECTS OF ADOPTING NEW PRECESSION AND EQUINOX CORRECTIONS ON THE FRAME 
OF THE CONVENTIONAL TERRESTRIAL REFERENCE SYSTEM 

1.1 Transformation Between Conventional Inertial (CIS) and Terrestrial 
Reference Frames (CTS) 

The transformation at an epoch T between the CIS at some fundamental epoch 
(e.g., 1950.0) and the CTS is 

[CTS] = SNP(M) [ CIS] (1) 

(see [Mueller 1981]). Here 

t 

s = M-* p ) M-y p ) Me) 
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is the earth rotation matrix, in which x p and y p are the polar motion components, 
and e is the Greenwich Apparent Sidereal Time (corresponding to the epoch T) 
computed from 

0 = (GMST)o + oj UT1 + Eq. E. 

6 

L. 

where (GMST)o is the Greenwich Mean Sidereal Time at 0 n UT1, is the conver- 
sion factor from mean time to mean sidereal time, and Eq. E. is the equation 
of the equinox (nutation in right ascension). The other matrices N, P, M in 
equation (1) are the nutation, precession, and proper motion matrices respec- 
tively [Mueller 1969, p. 123]. Parentheses around the M matrix indicate 
that proper motion is applied only in the case of a stellar CIS. 

Let prime (') denote the case with the precession, nutation and equinox 
changes introduced. The transformation equation (1) also holds for the cor- 
rected case: 

[CTS] ' = S'N'P'(M') [CIS] 1 (1‘) 

In this section only the precession and equinox changes are considered so that 
N' = N. From the definitions (or stipulations), one can determine directly or 
indirectly the relations between P' and P, M' and M, and [CIS] 1 and [ CIS ] at 
some epoch, leaving S' and [CTS] 1 to be solved for. 

One cannot solve for both S' and [CTS ] 1 simultaneously, hence some addi- 
tional constraint is needed. There are several options for the constraint, and 
they will be discussed later in Section 2.2. For the time being we will conform 
with the IAU adopted constraint, namely: Let the new ERP's be the same as the 
old ones at some epoch T u (in this paper T denotes the epoch, and t the time 
interval between T and some fundamental epoch, e.g. , 1950.0); solve for [CTS ] 1 
at this time, then keep it time invariant and solve for the resulting time 
variations in the new ERP's. 

1.2 The Effect in the Case of a Stellar CIS 

The new (1976) corrections for 1 uni solar precession in longitude and 
planetary precession in right ascension are [Williams and Melbourne 1981] 

APi - l'.'l/cy 

Ax - -0'.'029/cy 
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The correction to the equinox is Eo + Et, where Eo = Ot'525 is the offset at 
1950.0, E = l'.'275/cy, and t is the time elapsed from 1950.0 [Fricke 1931]. 


The new precession matrix P' can be written with sufficient approximation 
as 

P' = R 2 (Ant) R 3 (-Amt) P (2) 

with 

An = Api sine 
Am = APi cose - Ax 


where e is the obliquity of the ecliptic, and An, Am are the general precession 
changes in declination and in right ascension. Due to the equinox correction, 
the equation for the Greenwich Mean Sidereal Time is to change(without terms of 
higher order) to [Aoki et al . 1982] 

(GMST)q = (GMST) q + E 0 + Et (3) 

For the stellar (i.e., classical optical) CIS the change caused by the equinox 
correction at the fundamental epoch 1950.0 is 

[CIS]' « R 3 (-Eo) [CIS] (4) 

The new proper motion matrix is 

M' = Rz(-Ant) R 3 [(Am - E)t] M (5) 


The proper motion components in right ascension and declination are 

(u ) ' = (u ) + E - Am - An sina tan6 

Ur Ot 

(y 6 ) ' = (p 6 ) - An cosa 

Substituting the above new values of P', M' , [ CIS] 1 and (GMST)^(i ,e. , eqs. (2) 
- (5)) into eq. (I 1 ), one gets 

[CTS] 1 = R 2 (-x') Ri(-yp) R 3 [(GMST) 0 + o> e UT1' + Eq. E.] R 3 (E 0 + Et) N • 

• R 2 (Ant) R 3 (-Amt)P R 2 (-Ant) R 3 [(Am - E)t] M R 3 (-E 0 ) [CIS ] 

Except for (GMST) 0 » all rotation angels are small; neglecting the second- 
order terms, approximately, 

[CTS]' = Rz(-x') Ri(-yp) R 3 [(GMST) 0 + co e UT1 ' + Eq. E.] NPM [CIS] (6) 
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(For the above given ap x and E values, neglecting the modulation of NP will 
cause an error of less than O'.'OOOl in t * 10 yr.) Combining the above equation 
with the mentioned constraints at epoch T u : x^ = x p , y^ * y p , and UT1' ■ UT1, 
one obtains 

[CTS ] 1 - [ CTS ] 


The well-known conclusion is that in the case of the stellar CIS, the CTS 
and ERP's are unaffected because changes in the proper motion compensate for the 
equinox and precession changes. This statement is naturally valid not only at 
the epoch T u but at any time before or after. 

1.3 The Effect in the Case of a Non-Stellar CIS 

For any non-stellar (e.g«, VLBI or LLR) CIS, the proper motion matrix is 
no longer taken into consideration; the P* and (GMST)^ are the same as in the 
stellar case (eq. (2) - (3)). The relationship between [ CIS] 1 and [CIS] depends 
on the particular CIS under consideration. Generally, 

[CIS] 1 = Ej [CIS] 

If the considered CIS is aligned with the dynamic equator and equinox, then 
Ej = I, where I is a unit matrix. 

If the non-stellar CIS is aligned with the stellar system equinox at some 
epoch To, then Ej will be a little complicated. At this time due to the equi- 
nox correction, 

[CIS]! = R 3 (-Eo - Eto) [CIS] T (7) 

(More exactly, a second-order term could be considered.) The precession effect 
on the CIS's for the time interval t 0 between the fundamental epoch 1950.0 
and the alignment epoch T 0 is 

[CIS] T » P(to) [CIS] 

I 0 

[CIS]' = P' (to) [CIS]' 

I o 

With equations (2) and (7) one gets at the fundamental epoch 

[CIS] ' = R 2 (-Ant 0 ) R 3 [(Am - E)t 0 ] R 3 (-E 0 ) [CIS] » E. j [CIS] (8) 

i . e . , 
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Ej * R 2 (-Ant 0 ) R 3 [(Am - £)to - Eo3 (9) 

The corresponding corrections In right ascension (Aa F ) and declination (a«$ f ) are 

C I t I 


Aa F 
1 1 

A6 F 

E I 


* Eo + (E - Am) to - An to si no tani 

* -Ant 0 cosa . 


Now, substituting P\ (GMST) ^ * and [CIS ] 1 (i.e., eq. (2), (3) and ( 8 )) 
into eq. (r)» 

[CIS]' = Ra(-x^) Rx(-y^) R 3 [(GMST ) 0 + E 0 + Et + o» e UT 1 * + Eq. E] N • 

• R 2 (Ant) R 3 (-Amt) P Ej [ CIS] (1") 

As stated before, the ERP's are continuous, that is, at the alignment epoch T u , 

Xp * x p , = y p , UT1* = UT1. Thus 

[CTTS]' = R 2 (-X p ) Ri(-Yp) R 3 [(GMST)o + « e UT1 +- Eq. E] R 3 {E 0 + Et u ) N • 

• R 2 (Ant u ) R 3 (-Amt u ) P Ej [CIS] = 

= SNP R 2 (Ant u ) R 3 C E o + (E - Am) t u ] Ej [CIS] (10) 

If the CIS is linked with the stellar system equinox at epoch T 0 , i.e., Ej is 
expressed by eq. (9), then 

[CJS]' ± SNP R 2 [An(t u -t 0 )] R 3 [ (E-Am) ( t u -t Q ) ] [CIS] 

4 S R 2 [An(t u -t 0 )] R 3 [(E-Am)(t u -t 0 )] NP [CIS] (10') 

As pointed out previously, the modulation of NP is negligible, but the modulation 
of R 3 (e), included in S, must be taken into consideration. 

Rale) RaCAn(t u -t 0 )J = {R a (e) R 2 [An(t u -t Q )] R 3 (-e)} R 3 (e) 

* Ri[An(t u -t Q ) sine] R 2 [An(t u -t Q ) cose] R 3 (e) 

Substituting this into equation (10'), 

[GTS ] 1 ± Ri[An(t u -t 0 ) Sine] R 2 [An(t u -t Q ) cose] R 3 [(E-Am)(t u -t Q )] SNP [CIS] 

Thus for the case of CIS alignment with the stellar system 
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mV * fU[An(t u -t 0 ) Sine] RzUn(t u -t 0 ) cose] R3[(E-Am)(t u -t 0 )][CTS] 


(ID 


If the CIS is aligned with the dynamic equinox , that is, Ej a I, then 
[CTS] ' - SNP R 2 (An t u ) R 3 [E 0 + (E-Am) t u ][CIS] 

Thus 


[CTS] 1 = Ri(An t u Sine) R 2 (An t u cose) R 3 [E 0 + (E-Am) t u ] [CTS] 


(12) 


If the alignment is made over some time period (say, five days or so) T u is the 
mean epoch of alignment, and the values sine and cose are the mean values 
within this time span and can be averaged to zero. In this case 

[CJS] 1 = R 3 [(E-Am)(t u -t 0 )][CIS] (11') 


for the CIS linked with the stellar system equinox, and 

[CTS] ' = R 3 [Eo + (E-Am) t u ][CTS] 


( 12 *) 


when aligned with the dynamic equinox. Thus the relation between the new and 
old CTS's is a small rotation around the third axis. Expressed in longitude 
(positive to the East), 

6X = X 1 -X s (Am-E)(t u -t 0 ) (11") 

for the CIS linked with the stellar system equinox, and 

<5X = X 1 - X = (Am-E) t u - E 0 (12") 


when aligned with the dynamic equinox. 


For a CIS linked with the stellar system, if t u = t Q , then 6X = 0; 
otherwise a shift in longitude is necessary. As for a CIS aligned with the 
dynamic equinox, the CTS longitude origin shift generally cannot be avoided. 


1.4 The Effect of the Time-Invariant CTS on the ERP's 

The new CTS' at the time of alignment T u can then be determined as out- 
lined in the previous section, i.e., in the stellar CIS case [ CTS] 1 = [CTS] , 
and in the non-stellar cases as given by eqs. (11), (11 * ) or (12), (12'). 
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The next step Is to keep the new CTS time Invariant and to find the resulting 
ERP's at any time other than T u . Substituting eq. (11') for the left side of 
eq. (1"), and eq. (9) In the right-hand side, after some derivation and 
neglecting second-order terms, one gets 

[CIS] 4 Ra(-X*) Ri(-yp) R3[(CiMST)„ + u e UTl' + Eq. E] • 

* R 2 [an(t-t 0 )] R3[(E-Am)(t-t u )] NP [CIS] 

Comparing this equation with eq. (1), | 

R 2 (-Xp) Ri(-yp) R 3 C ( GMST ) 0 + u> e UT1 * + Eq. E] R 2 [An(t-t 0 )] • 

• R 3 C(E-Am)(t-t u )] = S 
or 

R 2 [-Xp + An(t-t Q )cos 0 ] RxC-Yp + An(t-t 0 )sin 6 ] R 3 {(GMST)o + Eq. E + 

+ [o) e UT1' + (E-Am)(t-t u )]} = 

= R 2 (-X ) Ri(-yJ RaC (GMST ) 0 + % UT1 + Eq. E] 

P P e i 

i 

From the above it is obvious that over a limited time span (otherwise second- 
order terms must be added), 

Ax p = x^ - x p = An(t-t 0 ) cose 

A y p = y p - y p = An(t-t Q ) sine (13) 

AUT1 = UT1' - UT1 = (Am-E;(t-t u )/w e 

The above are in the case of a non-steVlar CIS linked with the stellar system. j 

For the dynamic equinox alignment, substitute eq. ( 12 ') for the left side of 
eq. (1") and let Ej = I . The results are 

AXp = Ant cose ; 

A y p - Ant sine (14) 

AUT1 = (Am-E)(t-t u )/co e 1 

For both cases AUT1 is the same; so is the rate of AUT1: 

■ 

= (Am-E)/w e = -0.157 ms/yr (15) 

I 
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In conclusion, In the case of a non* stellar CIS, changes in the preces- 
Monti constant and the equinox will result In changes in both the CTS and the 
ERR 1 s » The CTS change 1 $ a longitude origin shift. The ERP changes are diurnal 
terms In the polar motion components with amplitudes linearly Increasing with 
time and a constant rate change In UTl, One point worth stressing Is that these 
are the differences of the same system (technique) between the new and old 
cases, not the differences between different systems (techniques). Also the 
diurnal tern which Is evident In polar motion Is not the diurnal true polar 
motion, but an artifact due to the time Invariant CTS constraint applied* 

2, GENERAL SOLUTION*. SMALL CIS ROTATIONS AND THEIR EFFECT ON THE CTS 

2.1 Changes in the Earth Rotation Parameters 

in the general case, eq. (13) or (14) can be written In the following 
form 

4 Xp * sine + eoso 

&y p * cose + «i sine ( 16 ) 

ae * «<s.j 

where the small angles represent the changes In the sense 

N'P'M' CMI' * Rita) M*s) NPM [CI$3 

Since 

9 * (GMST)q + w q UTl + Eq. E 

there are several possibilities for changing 9, If the nutation Is assumed 
to be unchanged, either may be absorbed Into (GMSTh, l.e,, It becomes a 
change In the Greenwich Mean Sidereal Time*, or, as before, It may go Into 
UTl (&UT1 * as MJi or It can be Incorporated partially In (GMST)q and par- 
tially in UTl, When Is placed (fully or partly) Into UTl, then If UTl 
Is still to be continuous at the epoch T u , the longitude \ has to absorb the 
one-time discontinuity as shown before, Finally, If as Is a nutation correc- ( 
fion, then 04 must be combined with Eq, E (see Section 3), 
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The corrections for precession, for the equinox, and for proper motion 
may be written in the following forms respectively 

ap * Rft(Ant) Rj(-Amt) 

Ej * ME Ijt ) RA£ h ) 

AM * R&(AM«) RjUMa) 

where, from comparisons with earlier results, in the case of the stellar CIS, 
AM* * "Ant, AM 3 * (Am-I)t, Ej * 0, E Ia » -E q j for the non-stellar CIS 
aligned with the dynamic equinox, E^ * E Ia * aM* « AM* « O', and in the case 
of the non-stellar CIS linked with the stellar system, E T * *Anto, E, » -Eg + 
(Am-E)to , AM* » AM * * Q« 


In any of the above cases 

«i * 0 

<x* * Ant f Ej 4 + AM* (17) 

« 3 • -Amt + eJ 3 + AM 3 

Thus, for example, In the case of the stellar CIS 

* «* * Q 

QL* •• -Eo ». Et (17') 

AO * Eo + Et 

If we let AO be the a(GMST)q» as we did before, then eq, (17 ') Is equivalent 
to eq, (3), 


(m 


In the case where the non-stellar CIS Is linked at To. with the stellar 

system, 

* 0 

a* * An(t-to) 
cig 81 -Eq + (Aro-t)ti - Amt 

i 

AO * Eq + (E-Am)to + Amt 

If we let 4(6MST)o a E 0 + Et, and let the ERR's be continuous at T u , then 
eq, (17") Is equivalent to eqs, (11") and (13). The analogy can also be 
established for the case of the non-stellar CIS linked to the dynamic equinox 
(eq, (14)), 
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2.2 Options to Change the CTS (Due to A0) 

As shown, in the case of equinox and precession constant corrections, 

A0 » Eo + Et for the stellar system, and 

A0 a -Ej a + Amt for the non-stellar systems 

A6 can also be written (assuming no change in nutation and the one-time 
discontinuity in UT1 absorbed in the longitudes, fix, mentioned earlier' 1 , 

A6 = a(GMST) 0 + co AUT1 + fiX 

c 

Thus, as stated above, one can absorb A6 either in a(GMST)o» or AUTI, or fix, 
or in some combinations of these. To get a definite (unique) solution, some 
constraint is needed. Mathematically, there are quite a number of possible 
choices for such a constraint, but practically only a few are meaningful. 
Below we deal with three sets of options. Which option is the best surely 
will be the subject of many discussions. 

Set A Options . Here the basic requirements are: (i) no discontinuity 

in ERP's at the epoch T u , (ii) the change in the Greenwich Mean Sidereal Time 
formula is the same for all CIS's, though different for each option. 


Set A Options 

Stellar CIS 
A0 = Eo + Et 

Non-stellar CIS 
A8 = -Ej 3 + Amt 

a(GMST) o 

Eo + Et 

Eo +^Et 

Option I co g aUTI 

0 

(4m-E)(t-t u ) 

5X 

0 

(Am-E)t u - E Ia - Eo 

A(GMST)o 

0 

f. 

0 

Option II w aUTI 

E(t-tu) 

Am(t-t u ) 

fiX 

Eo + Et u 

- e I 3 + Amt u 

a(gmst) 0 

Amt 

Amt 

Option III o> e aUTI 

(E-4m)(t-t u ) 

0 

fix 

Eo + (E-4m)t u 

• E Io 


Options I, II, and III above are similar to Tables 1, 2, and 3 respec- 
tively in [Williams and Melbourne 1981]. (The main difference appears to be 
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that for the non-stellar cases the general precession in right ascension Am 
is replaced by what they call "the average value over all observations of 
the effects of the precession corrections in right ascension" <a p >. For 
VLB I <&p> = Api cose - A& = Am, but for LLR <ot p > = Ap x - Ax t Am.) We 
already elaborated on Option I in Section 1.1. Option 1 is the presently 
accepted approach for the new FK5 CIS. But, as pointed out by [Williams 
and Melbourne 1981], for future possible new improvement of the precession 
constant and equinox corrections, this option might not be the best. They 
favor Option III because the space techniques are becoming the dominant source 
of Information about the transformation parameters between the CIS and CTS 
frames and because this option keeps UT1 invariant to improved values of the 
precession constant and the equinox position for the space techniques. The 
common geodetic disadvantage of Set A options is the required shift in the 
longitude origin (except in Option I for the stellar CIS case), the worst 
thing being that these shifts are different in the cases of stellar and non- 
stellar CIS's. 

Set B Options . Here the basic requirements are: (1) no change in the 
CTS, i.e., (5A = 0, (ii) as before, the Greenwich Mean Sidereal Time formula 
change is the same in all CIS cases, but different for each option. 

The major inconvenience of Set B options is the change in UT1, not only 
in the rate, but also in the necessary discontinuity. The value of the dis- 
continuity would need to be added with opposite sign to the UT1 at the epoch 
when the changes (new constants) are introduced. 
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Set B Options 

Stellar CIS 
A9 = Eo + Et 

Non-stellar CIS 
A9 * -Ej 3 + Amt 

sx 

0 

0 

Option IV a(GMST) o 

0 

0 

u> AUT1 

c 

E 0 + Et 

-E t + Amt 

1 3 

6X 

0 

0 

Option V A ( GMST ) o 

Eo + £t 

Eo + Et 

w e AUT1 

0 

-Er - Eo + (Am-E)t 

i 3 

6X 

0 

0 

Option VI A (GMST) g 

-E T + Amt 

1 3 # 

-E t + Amt 

* 3 

w e AUT1 

Eo + E t + (E-Am)t 

A 3 

0 


Set C Option . Here the basic requirements ares (i) no change in CTS, 
( i i ) no change in UT1, i.e., as is entirely absorbed in a(GMST)o. 


Set C Option 

Stellar CIS 
A9 = E 0 + Et 

Non-stellar CIS 
A0 = -E t + Amt 

13 

sx 

0 

0 

Option VII oj aUTI 

c 

0 

• 

0 

A (GMST) o 

Eo + Et 

-E t + Amt 

I 3 


Although this option is probably the preference of geodesists, it may seem 
to be unorthodox from the traditional astronomical point of view. How can 
the formulae for Greenwich Mean Sidereal Time for different CIS's be differ- 
ent? What will the astronomical meaning of (GMST)o be? However, one can 
view the formula for Greenwich Mean Sidereal Time as composed of two parts: 
The first part, (GMST) 0 , has its original astronomical meaning, while the 
second part, a(GMST) 0 is only a correction particular for a given CIS. It 
would make sense that since the changes in precession and the equinox affect 
different CIS's in different ways, this correction should also be different. 
From this point of view, Option VII seems plausible and even preferable for 
geodetic use. 
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It should also be noted that after the new equinox and precession 
changes are Introduced (once) into a(GMST) 0 , this option could become the 
equivalent of referencing the GMST to a fixed point on the equator, instead 
of to the mean equinox of date, the current practice. As pointed out by a 
number of authors, the advantage of such a change would be overwhelming and 
would make the future CTS stable against changes in the precession constant, 
etc. [Guinot 1979, Murray 1979, Williams and Melbourne 1981, Mueller 1981]. 

3. EFFECT OF ASTRONOMICAL NUTATION CHANGES ON EARTH ROTATION PARAMETERS 

According to the principle in Section 2.1, it is also easy to deal with 
any future changes in nutation. The nutation matrix N is [Mueller 1969] 

N = Ri(-e - Ae) R3 ( —Aip ) Ri(e) 
i Ri(-Ae) R 2 (Aip sine) R 3 ( -A41 cose) 

where Aip and- Ae are the nutation in longitude and obliquity respectively, and 
e is the obliquity of the ecliptic. If 6Ae and <$Aip are the respective cor- 
rections to Ae and Aip , then one can easily obtain the nutation correction 
matrix, 

AN - Ri(-6Ae) Ra ( <S Aip sine) R 3 ( -<S Ai|i cose) 

Thus in the notation of eq. (16), 

-6 Ae = cti 

6Aip sine - a 2 
-6Aip COSe = a 3 

and, therefore, 

AXp = 6Ae sine + 6 Alp sine cose 
Ayp = -<5Ae cose + 6 Alp sine sine 

A0 - -a 3 = 6 Aip COSe 

Thus, as expected, the effects on polar motion components are diurnal terms 
(5 Aip and 6Ae are long periodic). Again, this is a diurnal artifact in polar 
motion due to the introduction of the new nutation and not diurnal true polar 
motion. 

As far as the term A6 = 6 Aip cose is concerned, if it is incorporated 
into the Eq. E, neither the longitude origin nor the UT1 will be affected. 
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APPENDIX 

EFFECTS OF DIFFERENCES BETWEEN VARIOUS CTS'S AND CIS'S ON EARTH ROTATION 
PARAMETERS AND THE DETERMINATION OF SUCH DIFFERENCES 

The two CIS's (and two CTS's) inherent in two different techniques (e.g., 
SLR and VLB I) are generally not exactly identical [Mueller 1981]. Suppose the 
relation between the two CIS's at any epoch is (common nutation (N) and 
precession (P) matrices are assumed to be used in both techniques) 

[CIS] 11 * M“i) *»(«») M«i) [CIS] 1 (A.l) 

Similarly, the relation between the two CTS's is 

tml 11 = MBi) Ra(Ba) M3 3 ) C CJS ] 1 (A. 2) 

where a., and ^ are small rotation angles about the axes "i". 

The transformation from CIS to CTS again is 

[CJS] 1 = S 1 N P [CIS] 1 (A. 3) 

and 

[CTS ] 11 = S 11 N P [CIS ] 11 (A. 4) 

Substituting eq. (A.l) for the last term of the right-hand side of eq. 

(A. 4), and eq. (A. 2) for the left-hand side, 

Ri(Bx) MB*) MBaJCCIS] 1 = S 11 N P Ri(ai) R 2 (a 2 ) R^MCCIS] 1 

After some reduction, neglecting second-order terms, 

[CTS ] 1 = Ri(-6i + ai cose + a 2 sine) R 2 (-0 2 - a* sine + a 2 cose) • 

• R 3 (-B 3 + a,) S 11 N P [CIS] 1 

Comparing the above equation with (A. 3) 

S 1 = Ri(-Pi + ai cose + ct 2 sine) R 2 (-e 2 - a x sine + a 2 cose) • 

• R 3 (-B 3 + ci 3 ) S 11 
Or 

-*y p * -(yp-yj 1 ) = -Bi + <*i cose + a 2 sine . 

-ax = -(xj - xj 1 ) = -3 2 - a x sine + a 2 cose : 

r r r 
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« a AllTl = wjUTl 1 - UT1 11 ) = -0a + ct 3 (A. 5) 

6 6 

Thus the CTS differences (B angles) cause biases in all earth rotation 
parameters. Because of the modulation of the earth's diurnal rotation, the 
effect of CIS differences (cti, a 2 ) on polar motion components are diurnal 
terms, while the effect of a 3 on UT1 is again a bias. 

The direct way to determine all the 0 angles is the method of station 
collocation, i.e., to position two different types of techniques at the same 
location. 


The "observation" equation is 


= vl 


II 

$1 


-Xi =- 

62 

63 

+ 


0 

-0s 

02 


03 

0 

-01 


-02 

01 

0 



x i 


x i 



+ c 

y i 


2 i 


z i 


+ ^i 


where xjj and x?* are the determined coordinates of the same collocated station 
i in the two CTS's, 6 ,. is the translation vector, and c is the scale difference. 
One must have at least three collocated stations if all seven unknowns are to 
be solved for. 


For connecting the CIS's, there are a few methods such as the use of 
space astrometry to connect the stellar CIS and the radio source CIS, or 
using differential VLBI (which, for example, was used when the Viking Mars 
Orbiters and a quasar were near eclipsing) to connect the planetary and radio 
source CIS's (see [Kovalevsky and Mueller 1981]). These are direct approaches. 
One indirect method is via station collocation, i.e., using the earth as an 
intermediate body (see [Kovalevsky 1980]): First by station collocation one 

determines the CTS difference (0 angles) as above, then through earth rotation 
parameter differences determined over the same time period one finds the CIS 
difference (a angles). Eq. (A. 5) is the basi's for connecting the two CIS's. 
More details on this subject may be found in [Mueller et al . 1982]. 

When considering the above method one should note that the diurnal 
polar motion difference terms in eq. (A. 5) will show up as long as there 
are differences between the two CIS's (i.e., o 3 and ct 2 exist). This may even 


17 


be the case in situavons when one (or both) of the techniques solve for 
rotations of its CIS, resulting in no (individual) diurnal polar motion. 
This, of course, would mean that the adopted precessional constant is 
discarded. 
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2*2 Utilization of Range-Difference Observations 
in Geodynamics (Research Contract NAS5-25888) 


2.21 Utilization of Simultaneous Lageos Range-Differences 
in Geodynamics 

Introduction 

The following is a summary of the research performed during 
the past six months under the lageos project, dealing with the 
utilization of simultaneous laser range-differences (SflD) for the 
determination of earth orientation and baseline variations. 

Reported are some results from the Aug. 1980 Lageos data 
collected during the short MERIT campaign, and simulations for a 
possible station arrangement for the main campaign (to begin in 
1983) . 

2.211 Simulations for a proposed MERIT83 laser network. 

Based on an optimal global laser station distribution (likely 
to be realizable by mid- 1983) proposed at a recent meeting of the i 
study group (cf . COTES proposal in last semiannual report) , a 
simulation study for baseline recovery was performed. Except for 
the fact that different stations (seventeen total) are involved, 
this simulation was similar to the one previously reported for 
the MERIT80 network in the last report. The station locations 
and the data distribution are given in Table 1 [ (a) -(b) ]. 

Baseline estimates and their statistics were computed for both 

4 

the range and the SRD adjustments. In order to assess the effect 
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of orbital biases on the baseline recovery, the orbit used in the 

adjustments (range and 3RD) was biased as follows : 

Radial bias : 2.00 m 

Along track bias : 0.60 m 

Across track bias : -1.20 m 

Two different adjustments were performed. In the first case 
the coordinates of all stations were obtained in a simultaneous 
adjustment based on the data collected from all baseline pairs. 

On the basis or this solution the baselines between all possible 
station combinations were obtained along with their formal 
accuracies and differences with respect to their "true” values. 
The results of this solution for the station coordinates are 
given in Table 2 for the range adjustment and in Table 3 for the 
SBD adjustment. The baseline results are shown in Table 4. 

As it can be seen from the last table, in all cases except 
for two, the baseline lengths have been overestimated although 
the errors in the SRD case are about an order of magnitude 
smaller than the ones for the range adjustment. Since the radial 
bias results in an ''expansion" of the network of satellite 
positions, this should come as no surprise. The stations have a 
global distribution and since the observations from all stations 
are adjusted simultaneously, their positions become 
interdependent and the aforementioned expansion affects all of 
them similarly. 

The results of this first adjustment prompted us to test the 
recovery of baselines from individual adjustments. In this 
second case the data collected from each pair of stations are 
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adjusted independently and the estimated baselines are only the 
ones defined by coobsecving station pairs. The results of this 
second type solution are shown in Table 5. What is obvious again 
is that the SRD results are again superior to the range results 
for baselines and station coordinates alike. The quality of the 
results with respect to the latter is characterized by the norm 
||X(| of the six coordinate differences between the true and 
estimated positions of the stations defining each baseline. 

The most interesting observation though in this solution is 
that on the basis of the same data, the range adjustment now 
underestimates the baselines and the recovery errors are all 
negative. For the SRD results, there seems to be no bias 
preference and those errors are rather randomly distributed and 
in almost all cases at the centimeter level. The three baselines 
for which the range adjustment has given better results than the 
SRD, all have lengths in excess of 7000 km and very few 
observations. As it has been previously reported the SRD mode is 
much more geometry dependent than the range mode, and as the 
results of this table show it admits of its limitations very 
eagerly (note the formal accuracies on those baselines !). 

Unlike the SRD mode, the formal accuracies for the range mode 
give no hint whatsoever as to the real accuracy of the results. 
Even though the recovery errors are of the order of a few 
decimeters in all cases, the reported <r*s are hardly ever higher 
than 2 cm ! 

On the basis of these simulations one can conclude that the 
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Table 2 Recovered Station Coordinates (Range Mode) 


STATION NO. J 7051 
APRIORl ESTIMATE 
ADJUSTMENTS 

X 

-2516274.096042 

-1.287283 

AOJUSTEO POSITION 

-2516276.183323 

STATION NO. t 7063 
APRIORl ESTIMATE 
ADJUSTMENTS 

X 

1130304.017676 

-0.267394 

AOJUSTEO POSITION 

1130304. 5502 fH 

STATION NO. : 7069 
APRIORl ESTIMATE 
ADJUSTMENTS 

X 

961533.600910 

-0.264610 

AOJUSTEO POSITION 

961533.336300 

STATION NO. i 7086 
APR10RI ESTIMATE 
ADJUSTMENTS 

X 

-1324510.442373 

-0.902349 

AOJUSTEO POSITION 

-1326511.344721 

STATION NO. I 7090 
APR IORI ESTIMATE 
ADJUSTMENTS 

X 

-2389125.331291 

-0.975459 

ADJUSTED position 

-2389126.306750 

STATION NO. » 7091 
APRIQRI ESTIMATE 
ADJUSTMENTS 

X 

1492212.741998 

-0.037315 

ACJUSTEO POSITION 

1492212.704682 

STATION NO. t 7095 
APRIORl ESTIMATE 

adjustments 

X 

3392750.871854 

1.469007 

ADJUSTED POSITION 

3392752.340862 

STATION NO. : 7120 
APRIORl ESTIMATE 
ADJUSTMENTS 

X 

-5464096.682969 

-2.035242 

AOJUSTEO POSITION 

-5464098.718211 

STATION NO. * 7901 
APRIORl ESTIMATE 
ADJUSTMENTS 

# 

X 

3844341.318863 

1.470786 

ADJUSTED POSIT ION 

3844342.789649 


V 

Z 

-4198843.469479 

4075154.588717 

-1.351438 

1*822092 

-4198844.820917 

4075156.410809 

Y 

Z 

-4831721.449137 

3993759.624496 

-1.937618 

1.854899 

-4831723.386754 

3993761.479395 

Y 

Z 

-5674186.967561 

2740519.740502 

-2.204506 

1.526215 

-5674189.172147 

2740521.266717 

Y 

Z 

-5332139.932091 

3231791.055906 

-1.830074 

1.691947 

-5332141.762966 

3231792.747852 

Y 

Z 

5042839.037557 

-3078750.720221 

1.621840 

-1.359841 

5042840.659397 

-3070752.088062 

Y 

Z 

-4458121.790935 

4296005.408571 

-1.836603 

1.862885 

-4458123.627538 

4296007.351456 

Y 

Z 

783278.256725 

5325906.606633 

-0.126506 

2.043665 

783270.130219 

5325900.650490 

Y 

l 

-2402363.153199 

2240358.272655 

-0.395526 

1.552528 

-2402363. 54G 725 

2240359.825183 

Y 

Z 

-134247.357044 

5070549.689834 

-0.208830 

1.992617 

-134247.565874 

5070551.682451 
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STATION NO* I 7907 
apriqri estinatb 

ADJUSTMENTS 

X 

1941330*114713 

-0.233477 

V 

-180*024 .122161 
- 2.273*46 

2 

-1774312.783770 

0.774573 

Aiijusrao position 

misao.aiuu 

-580*026.373407 

-1776312.211177 

STATION NO. l 77U 
APMQIU BSTIMATB 
AQJUSTMBNTS 

X 

mmuyum 

U37063* 

V 

0.0 

-0.343076 

l 

49339S0. 633358 
2.103731 

ADJUSTED POSITION 

4Q22037, 198262 

-0.343076 

4935852.737270 

STATION NQ» \ 7714 
APRIQRI BSTINATB 
AQJUSTMBNTS 

X 

4074413.IQ48T7 

U84T344 

V 

731763.678222 

-0.187477 

* 

4801472,271034 

2.037665 

ADJUSTED POSITION 


731763**88523 

4801474.330677 

STATION NO. X 7m 
APRIQRI ESTIMATE 
ADJUSTMENTS 

X 

-AiamT.'mSfiT 

- 1.710228 

Y 

3*20176.370484 

1.738147 

1 

3637871,317704 

1.771707 

ADJUSTED POSITION 

-41*1437.507815 

32*0178. 108633 

3637873.111614 

STATION NO. I 7740 
APRIQRI 6STIHATB 
ADJUSTMENTS 

X 

4728437.290679 

1.733163 

Y 

1710493,461733 

-Q.08QS64 

* 

3817377.771472 

1.757693 

A0JUST8Q POSITION 

4728638.783943 

1710473,381171 

3817379,751187 

STATION NO. X 7742 
APRIQRI ESTIMATE 
ADJUSTMENTS 

X 

4350757.298444 
U 6 166 10 

y 

437367,304711 

-0.203136 

2 

4408076.773477 

2.0*2975 

ADJUSTED POSITION 

4350760.875084 

637347,301373 

4408098, 7764T4 

STATION NO* X 7743 
APRIQRI ESTIMATE 
ADJUSTMENTS 

X 

-4843816, 65328T 
-0.801772 

1 

1343390,081748 

2.071760 

* 

-4488060.975096 

-1.331585 

ADJUSTED POSITION 

**4263817.433037 

1943332.733708 

-4488062.306640 

STATION NO. 1 7999 
APRIQRI ESTIMATE 
ADJUSTMENTS 

X 

4 13003 U 4878 74 
1.644677 

y 

1106638.602427 

-0,043308 

* 

4714882. OT47S0 
1.778227 

ADJUSTED POSITION 

4130033.134372 

1106638,538717 

4714884,073187 
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Table 3 Recovered Station Coordinates (SRD Mode) 


STATION NO. I 7051 
APR 10RI ESTIMATE 
ADJUSTMENTS 

X 

-2516274. 896042 
-0,562246 

V 

-4193843.469479 

0.077679 

2 

4075134.588717 

0,532835 

ADJUSTED POSITION 

“2516275.458287 

“4198843.391300 

4073135.141532 

STATION NO. J 7063 
APRI0R1 ESTIMATE 
ADJUSTMENTS 

X 

1130304.817676 
=>0, 412984 

Y 

“4831721.449137 

-0.339656 

2 

3993739.624496 

0.568841 

AOJUSTEO POSITION 

1130104.404692 

-413 1721V 780792 

3993760.193337 

STATION NO. t 7069 
APRIORI ESTIMATE 
ADJUSTMENTS 

X 

961533.600910 

-0.493799 

Y 

-5674186.967561 

-0*372649 

2 

2740519.740302 

0,304026 

AOJUSTEO POSITION 

961533.107111 

-5674187,340210 

2740520.244528 

STATION NO. t 7086 
APRIORI ESTIMATE 
ADJUSTMENTS 

X 

-1324510.442373 
**0. 614059 

Y 

-5332139.932091 

-0.110826 

2 

3231791.033906 

0.312840 

ADJUSTED POSITION 

-1324511.056431 

-5332140.042917 

3231791.368746 

STATION NO. s 7090 
APRIORI ESTIMATE 
ADJUSTMENTS 

X 

-2389125.331291 

0,207588 

Y 

5042839.037557 

0.480930 

£ 

-3078730.728221 

0,222993 

ADJUSTED POSITION 

-2389125.043703 

5042839.518487 

-3078730.505229 

STATION NO. i 7091 
APRIORI ESTIMATE 
ADJUSTMENTS 

X 

1492212.741998 

-0.348903 

Y 

-4458121.790935 
- 0.347854 

2 

. 4296005.480571 

0.591585 

ADJUSTED POSITION 

1492212.393095 

-4458122.138789 

4296006,080137 

STATION NO. » 7095 
APRIORI ESTIMATE 
ADJUSTMENTS 

X 

3392750.871854 

0.250456 

Y 

783278.256725 

-0.109695 

2 

3325906.606633 

0.704342 

AOJUSTEO POSITION 

5392751.122310 

783277.947029 

3323907.310973 

STATION NO. i 7120 
APRIORI ESTIMATE 
ADJUSTMENTS 

X 

“5464096.682969 

“0.589684 

Y 

-2402363.153199 

0.486807 

2 

2240358.272655 

0.474197 

AOJUSTEO POSITION 

“5464097.272654 

“2402362.666392 

2240358.746033 

STATION NO. « 7901 
APRIORI ESTIMATE 
ADJUSTMENTS 

• ■ 

X 

3844341.318863 

0.175767 

! Y 

“134247.357044 
*>0.397336 

2 

3070549,689034 

0.686721 

AOJUSTEO POSITION 

3844341.494631 

“134247.734380 

5070530.376333 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


STATION NO. I 7907 
APR l OR I ESTIMATE 
ADJUSTMENTS 

AOJUSTEO POSITION 


STATION NO. X 7911 
APRIORI ESTIMATE 
ADJUSTMENTS 

ADJUSTED POSITION 


STATION NO. I 7914 
APR IQRI ESTIMATE 
ADJUSTMENTS 

ADJUSTED POSITION 


STATION NO. x 7933 
APRIQRI ESTIMATE 
ADJUSTMENTS 

ADJUSTED POSITION 


STATION NO. X 7940 
APRIQRI ESTIMATE 
ADJUSTMENTS 

ADJUSTED POSITION 


STATION NO. X 7942 
APRIQRI ESTIMATE 
ADJUSTMENTS 

ADJUSTED POSITION 


STATION NO. J 7943 
APRIQRI ESTIMATE 
ADJUSTMENTS 

AOJUSTEO POSITION 


STATION NO. a 7999 
APRIORI ESTIMATE 
ADJUSTMENTS 

ADJUSTED POSITION 


X 

Y 

2 

1941330. U4913 

-5802024.122161 

-1796312.985770 

-0.444374 

-0*522985 

0.256872 

1941329.670339 

-5802024.645146 

-1796312,728899 

X 

Y 

2 

4022035.767630 

0.0 

4933530.635358 

0.198181 

-0.407657 

0,673034 

4022038*965012 

-0*407657 

4933551.308392 

X 

Y 

l 

4074613.304579 

931963*678222 

4801492.271034 

0.301324 

-0.373080 

0.668115 

4074613.606103 

931963.305142 

4801492.939149 

X 

Y 

2 

-4121637,799587 

3220176.370484 

3637871.319704 

0.044690 

0.500912 

0.490394 

-4121637.754897 

3220176,931396 

3637871.810098 

X 

Y 

2 

4728637.250678 

1910493.461735 

3817397,791492 

0.444574 

-0.387945 

0.619332 

4728637.693252 

1910493.073790 

3817398.410824 

X 

Y 

2 

4550759.250444 

639567.504711 

4408096.973499 

0.296323 

-0.433382 

0.648499 

4550759.554768 

639567.071329 

4400097.621998 

X 

Y 

2 

-4245816.653287 

1545350.031948 

-4488060.975056 

-0.249451 

0.494039 

0,067337 

-4245816.902774 

1545351.373987 

-4488060.907719 

X 

Y 

2 

4130031.489874 

1106630.602427 

4716882.074950 

0.324886 

-0.372003 

0.668722 

4130031.814759 

1106638.230423 

4716882.743680 
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SRD mode will in All likelihood provide more meaningful results 
in the presence of unmodeled orbital Jbiases than the range mode, 
and it will also give more reliable accuracy estimates for those 
results. Comparing the batch (global) solution to that of 
individual adjustments, the latter seems to be by far a better 
approach in the case of SRD observations, although the opposite 
is true for the range observations, compare for instance the 
level of recovery errors between Tables 4 and 5. 

2.212 Preliminary Results from Iageos Data Analysis. 

Lageos ranging data collected from ten stations over the 
period August 14-29,1980 (during the short HfiRIT campaign) were 
used in GE0SPP8 1 for baseline recovery. A total of 24240 ranges 
were selected with effort to balance the distribution among 
stations whose observability performance shows wiid variations 
(cf. station 7050 with over 60000 ranges during August, and 
station 7092 with hardly over 3000 in the same period of time) . 
The summary of the data distribution per pass per station is 
given in Table 6 [ (a)-(j) ]. The ill-conditioning of the normal 
equations due to the lack of origin of longitudes definition is 
overcome by applying a small weight in all three coordinates of 
ail stations, corresponding to a 0'= + 50 m. This way the origin 
of longitudes does not depend on a single station but rather the 
ensemble of them. The separation of the X and Y coordinates is 
thus not as good as it would be if one longitude were fixed 
absolutely, but that has no effect on the baselines. This high 
correlation between X and Y is also reflected in the estimated 
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formal accuracies for these coordinates/ Table 7. The orbital 
model and the constants used in the solution are shown in 
Table 8. Baseline results of the adjustment are given in 
Table 9 f and an analytical breakdown of the residuals after the 
adjustment are given in Table 10 [ (a)-(j) ]. Notice that the fact 
that station 7907 (ABELAS) is the one with the fewest 
observations (only 489) shows very clearly in the estimation of 
baselines which emanate from that station (Table 9). 

Care should be taken i,n comparing these results with other 
solutions for the fact that these baselines are reckoned between 
the optical centers of the corresponding laser instruments and 
not the stations* validation points. 

This investigation is now being completed, and the final report is in 
preparation by E. Pavlis, to appear in the report series of the Department 
of Geodetic Science and Surveying, The Ohio State University. 
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OBSERVATION SUNMAKV 

• mMUIUUMMII»UH>U«ilU»aiHI 


STATION IDENTIFICATION NU. I TO<*3 


NUMBER OF PASSES TRACKED I 


10 


NUMBER OF OBSERVATIONS « <i2BN 


'ft 


i 

* 

» 

9 

10 


ftEGINNlftG OA Jfc 

yyhhoo ****&£•* 


ENDING DA ft 
YYHMOQ HHHHSS.S 


aouaiA 

taut 

600427 
a 068?5 

ssasfe 


rood, 


4.Q 
614 6*o 
90040*0 

122631.0 

62146.0 
Idaho, u 

90332.0 

160010.0 

Vf«f?:8 


588552 

6 00422 


600814 16323O.0 

' 84127.0 

92334.0 
123440.0 
901 6.0 

...... Hli 42:8 

600627 16251 i.O 
600627 200656.0 
0 00626 62441.0 


8 00625 
6U082: 
60062 


DURATION 
SEC UNQS 


2410.0 

*185:8 

2444.0 


Q6SERVATiC«aS 


(UNE PUlN 


' UU A , t . 

PfcK A StCS* 


477 
202 
. 6 
639 

1350 

u« 


323.00 

3.44 

7.40 

261.50 

2.73 

0.0 

1.61 

376.76 

162.14 

*•1# 


(b) 


OBSERVATION 

■MUMMiaiM.Maaiiaaaa 


SO M AK I 

mmmmmmmmmmMMMu 


STATION luENTlFIGAT JUN NU. 

- ■ 1 ■» — '■ " i. ■ i ■ ■■■ — 


70 90 


N4JH4ER OF PASSES TRACKED : 


30 


NUMtiER Of UdSeRVA T1UNS 1 4143 


PASS 

NO* 


4 

5 

$ 

6 
9 

10 

11 

12 

p 

i 

19 

20 

& 

3 

25 

it 

26 

29 

30 


(Beginning date 

YYHHOO jhwNHSS.S 


ENDING DA Tt: 

Y YHHOO HHHHSS.S 


0UKAT1UN 

SECONDS 


OBSERVATIONS 


DENSITY LAO A 
(UNO POINT PtK A StCSA 


400614 

#00614 


73524.0 
UO 144.0 


400614 >1133 6.0 
400814 210657.0 


Milt 

600615 

600616 

600614 

600616 

400419 


600621 


93631.0 

162235.0 
1941 3.0 

91 ‘***3.0 

124 14 1.0 
i9Uo3*.U 

111332.0 

175013.0 
4.2*21.0 

629 6.0 

200 224.0 

63132.0 

i25&?:8 

171326.0 
*i»4 73o.O 

64746.0 

122010.0 

1334 1.0 
1647 3.0 
1U3o52.0 

143042.0 
172/46.0 

210429.0 

161 247.0 
1937 5.0 


800614 
600614 
600614 
400614 
600813 
600615 
d 00615 
6 00616 
400616 
600616 
*00819 
60O619 

HUUttiV 

600420 

400620 

600821 

500822 

600622 

600622 

600622 

600626 

800626 

600626 

600626 

400627 

600427 

*00427 

400427 

*00426 

600626 


6 1120 .0 
1136 6.0 
1617 7.0 
2148 tt.O 
102222.0 

165245.0 

202636.0 
9*642.0 

1259 1.0 
1953 1.0 

114926.0 

143136.0 

nmi: s 

204123.0 

91337.0 

74943.0 

111636.0 
17565 7.0 
2126 5.0 

920 7.0 

123626.0 

155759.0 

193427.0 

112919.0 

143215.0 

181144.0 

213936.0 

164735.0 
20 22 5.0 


2156.0 

4*2.0 

2519.0 

2471.0 

2631.0 

1110.0 

2653.0 

2037.0 

1040.0 

2767.0 

2036.0 

2461.0 
574.0 

1267.0 

2339.0 
2o*5 .0 

2491.0 

2409.0 

2421.0 

1090.0 

1436.0 

2642 .0 

1947.0 
93.0 

2636 .0 

2109.0 

2066.0 
2700.0 


iff 

Iff 

196 

141 

26i 

11? 

1?I 

136 

2U3 

18 

136 

104 

162 
173 
136 
155 
4 1 
66 
233 
115 
6 

169 

134 

214 


22.23 

13*u? 

13*6* 

11*94 

13.42 
12 . 6 * 
Uia3 
17.12 
13.3^ 

16.30 
14.94 
12*^2 
11. 40 
44.36 

17.20 

23.43 

43.96 
14*Uo 
1*.4U 
17.71 

1D.62 

26.39 

16.34 

12.20 

16.93 

13.30 

13.96 
13.o9 


(c) 


□ 8 SERVATION S U M M A R V 

sacBiiKaixiiKxiixKCBRiBiBaiiKiaR BBaaBsaac 


STATION IDENTIFICATION NO. ; 7091 
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YYHHOO ttlrtHSS.S 
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Y YHHOO HHHHSS.S 
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OBSERVATIONS 


density lag X 
<u tit point fek a secsi 


s jssunmftf 

" 588:15 l Hff?:8 Httff 

IH 

6.36 

*•#6 

l 58858 &SM 

888511* ilomio t*»’:8 

m 

7.63 

4.46 

(d) 

OBSERVATION SUMMARY 

■ i ■■■■■■■■■ « ■ i ■■■■■■■■■■■■■ ■■■Tfari mi 



STATION IuENTIFICaT Lki NU. S 

A) 92 NUMBER OF PASSES TRACKED t S 

NUHdfcK 

OF UBSEHVATIUNS l a2SA 


PASS 

NU. 


dcGINNIMS DATE 
Y YHHOO tnHHSS.S 


ENU1NG DATE 
Y YHHOO HHHHSS.S 


dura i ion 
secunos 


U6SEKVATIUNS 


DENSITY LAG A 
(ONE POINT PCK A SCLSI 


600623 133228.0 
600624 141333.0 
400624 173125.0 
600625 1613 4.0 
600#26 150956.0 


600*23 160149.0 
800624 1433 6.0 
600624 18 11 5.0 
6 00625 165055.0 
600626 151250.0 


1761.0 

1133.0 

2360.0 
2271 .0 

17* .0 


322 

266 

1273 

363 

9 


5**7 

*.03 

1**7 

o.2o 

19.33 
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STAtlUN 10ENTIPICAT KM NU. I 
* 1 — " * %* " ' M * 


UI^HVAUON SUMMARY 

naaa« «««»«» rim *««aa.B.«aia*aMa« ■*«****«*. 

7096 NUHdeR UF PASS** fHAClUU S 


6£UXNNlM# DA It 
VYMHUO *V4H**.S 


NUM48K UR U4$fckVAfiUNS 4 


8 N 01 NG OAT* 
Y YMMOO NHMHSS.S 

imp tin , J:s 

S 88 Si $ 

S 88 SH ittSiiS 


l &:8 

ij&u 


u 4 $i*?ArioN* 


otNAi rv 

iUNt puini 


LAW * * h 
PtK A ScC*t 


OUUVniDN MM A»y 

BaaaaaaaaaaaBaaaaaaaaaMaaaaaBBBBAaaaaaaa 


STATION IOfcNTlPlCAI ION NO* » 7114 


NUHttEK UP PASSES IHACKtO t 


NUMtftA UP URSIcNVAUONS I loo* 


•EMNNIftt DATE 
VVNMOU (WHHSS.S 


!S“i! 

*00421 

*00421 

ml 

•00426 

* 00*21 

aa 


il%4 6.0 
*210 3.0 
102434.0 
2u>140.0 

4&?:8 

BJ}«:8 

AJ07S4.0 


£ NO INS DATE 

VVNNDO HHHHSS.S 


sum mi m 

SWAZI II 06*. 9.0 


kuu«<ti iUW/.o 

sum mm 

SUSS it * 4 iiii :8 
SSSii ? | 88 fs 8:8 

*00427 211*14.0 


tlffcSS 

2101 .0 
. *47.0 
10** .0 
*76.0 
1040.0 
400.0 


iMSfckVAUUNS 


jUNt^PUlNI RtK U A *fcW 6 | 

6.40 ~ 

*6.62 

i4::5 

4*4i 

ul:i3 

o*r># 

1*4 .2V 

IOO.WO 


N3/ 04 M ft V A T 10(4 SUMMARY 

STATION IDfcNriFlCAnON NO. I 7115 NUKB6R OR PASStS TRACKtl) 3 U 


NUM48R Of UttSfcRVAUUN* S .2b* 


44U1NN1 NO OATc 
YYHMUU lr**MMSS. S 


SB 

* 00 * 1 * 

400*14 

*00414 

400414 

mi 

• 00*21 

400422 


111514.0 
1**114.0 
1141 1.0 
2l**4*.U 

110621.0 

143«**.0 

211024.0 
9* 7*4.0 

111219.0 
11462 7.0 

140436.0 
1244 7.0 


bNOlNO OATt 
YYMHOQ HNMMSS.S 


d 00614 
40041* 
400*15 
600*14 
•0041* 
400416 
400419 
400«19 
400*20 
• 00*21 
600422 


133*42.0 
200*29.0 
12 1624.0 

221246.0 
11%414.0 

150044.0 

214154.0 

101244.0 
1144*2.0 
1233*4.0 

143124.0 
12*4*9.0 


DORA f ION 
56CDNU6 

JHW " 
2001.0 
1021.0 

1346.0 

449.0 

? i »:8 


04S8RVAUUNS 


D8NSUY LAW X 

ION* ruin 7 Pt* a stu*) 


RMRR.RaLanMlM.MU.aMRl aa 2 a£«£a£ Y 


STATION XOtNUFluAr 4JN NO. 3 7120 


NUMbCR OF PASStS IRACAfcD s 


NuMdfcK Uir UoStRVAllONS * I9u% 


6EG1NN1N3 OATt 
YYMMOO ttmrtSS.S 


*00*14 

400*14 

*0041* 

40041* 

•0041* 

i&lil 

asa 

•00426 


131936.0 
46*434.0 

2*247.0 

121219.0 
l*223o.O 
1414 * 6.0 
12**43.0 

163425.0 

ssn:8 


8 NO 1 NO (MTfc 

YYMMOO MHMMSS.S 


•00414 

400414 

• 00*15 

• 00615 
400*1* 
• 00*21 
4 00422 
400422 


• 00*23 
400426 


135240.0 

171227.0 
3142%. 0 

122237.0 

160610.0 
1*0035.0 
133 636.0 
16*250.0 

mm 


00* AT ion 

SfctONU S 

1964.0 

429.0 
1297 .0 

614.0 

ftiiS 

iSli.U 

•65.0 


OSSENVaUUNS 


Iu ^ PuIni PtK^A *£Ci» 
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(0 


UlSERVATIQh 

mtmmmmmmmmmmmmmmmmmmmmmmm 


i U H H ARY 


STATION IDENTIFICATION NU. t 7907 

1 1 1 " ' 1 1 1 


N0M8ER OF HASSES TRACKED t 


20 


NUN1ER UF URSERVaUUNS * *o9 


PASS 

MJ* 


•eoinnin* date 

YYMHOQ rvtNNSS.S 


END INC DATE 

VVNNOO NHHMSS.S 


DORA HON 
>ECUNoS 


OBSERVATIONS 




k 

i 

* 

5 

* 

i 3 

0 

u 


ttttts 8^:4 

»00»t> - 

isii! 

88814 

•00119 
•00620 

Nil 



8 . 

BS&t 

234237* 7 
92152*3 
75452*3 
633 7.3 
101022*6 
1*432.3 
71730.1 
60230.0 

* 85:4 

vsotglu 

88K 


•OOll4 
lOOli* 
•00613 
• 00616 

sssst? 

i 00616 
• 00619 
i 00619 
600*20 

822811 

222211 

22221 s 

.00.2} 

.00.27 

. 00 . 2 . 

. 00 . 2 . 


Jtfji 

NBW:i 

• 0332*6 
6 *.313. 1 
102945*2 

m»!& 

60030*1 

ll&ki 

815 y i:t 

«a 

•»*22.* 


4 ! 5 ta 

io21*« 

360*2 

•92*9 

1162*6 

916.9 

961*6 

360.1 

x ii8:8 
620.0 
669 *0 

iStS:$ 

90.0 

?«0.0 


12 

St 

19 

32* 

8 

n 

9 

H 

29 


ll.Vl 

36*09 

n:W 

flsU 

iiiii 

*8:3S 

66.26 

60 . 3 V 

62.60 

61.6? 

36.66 

36.16 

6>.oa 

26.90 


(j) 


naccou.Ti n 


c tj m m * u y 


iMMUMaaiiwiMaMiuaiuiuiiKiuMUu 

STATION IDENTIFICATION NU* S 7963 NUMBER OF PASSES (RACKED t 32 


NUMBER UF UlStKVAUON* 1 3* lo 


FASS 

NU* 


•EOiNNiMt DATE 
YYNHOU TtiHHSS.S 

^ws* 

126959.9 
lt>03 7.7 
•1963 7*6 
*1344.9 
114022*1 
16652 2*7 
101513*0 
10io52*tf 
136030*6 
1701X3*0 
90622*7 
1217 0.2 
133013*1 
191022*3 
110013*2 
123030*2 
1507 7.7 
1*3030.1 
Ul2 622.7 
136232*6 
170913*1 
92930.2 

66132*0 
113322*7 
130922.7 
92537*5 
1236 0*0 
13572 2*7 


ENUlNO DATE fc 
VYMMOD HHHMS5.S 


DURATION 

SECONDS 


OBSERVATIONS 


UtNSin LA 0 X 
(ONE PUlNI HER A Seoli 


i 

3 

6 

5 

6 
7 
• 
9 

10 

11 

81 

15 

i 

* 


30 

31 


100615 

£<S«615 


?<*Hv815 
-vj.tfO* 16 

•00616 


•110616 

•09617 

•00617 


•00617 
•00616 
•00616 
•006 1 6 
lOOolo 
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Table 9 


6AS8L1N& bStlMAfeS KcLAlfcO STATISTICS > 


fcASfcLlNfctf 

STaIION* bTAlluN* 

APKlOHl fiST. 

A0JUST80 VAL. 

OlPF.IA-sJ 

51UHA 

KccAt IVt ACC 

k 

7 0*3 

• 4 

7u90 

U645951.161 

12646950. U47 

-o.vu 

0.0144 

3.3*U-0V 

2 

7043 

«*> 

7091 

sOkO^.lAi 

602032.169 

0.0** 

U.Oio 

1 .**0-0? 

3 

7043 

**> 

7094 

1000J 296.515 

10003295.633 

-0.3*2 

0.025 

6 *030-07 

6 

7043 


7099 

9896473.055 

9896471.52b 

-1.32* 

0.022 

3.1VO-OV 

» 

7043 

«•> 

7114 

3 562138.713 

356*137 ,642 


0.0*1 

2. 1*0-0* 

4 

7043 

*4 

7113 

isoisva.na 

3501891.797 

-1.3*1 

0.007 

2.320-0* 

7 

7043 

mw> 

7140 

7*44020.742 

7244019.261 

H.**2 

0.02* 

V.230-OV 

4 

7043 

*9 

7907 

5 9< 8036. 951 

5928019 <003 

-17.9** 

0.0*3 

3. *10-0* 

9 

tik>3 

»4 

7943 

U 106539.054 

1210653*) < 064 

-0.790 

O.olb 

3 .**0-0V 

10 

7090 

«"> 

7091 

12aJ01*U.062 

12638161)4,219 

0.13* 

0.02* 

* .320-OV 

Ik 

7090 

•4 

7o92 

6674009.770 

6674008.743 

-1.021 

0.02* 

*•*10-09 

u 

7090 

»4 

7099 

7*47620.43* 

7247520.743 

cull 

0.023 

7.5*0-09 

13 

7090 

•4 

7119 

U 768618.0 14 

117666U.337 

0.323 

O.Ul* 

3. *70-09 

14 

7090 

■ 4 

7113 

11810628.856 

11810629 .OU 

O.bH 

O.Olto 

3.200-09 

13 

7090 

■4 

7140 

9656456.579 

965643U.910 

9.331 

0.u21 

5.120-09 

14 

7090 

■ 4 

7907 

U 750456. 119 

11750658.620 

2.300 

0.03* 

3.920-09 

17 

7090 

>4 

799 3 

3196328.733 

3196328.646 

-o.o*? 

0.0*1 

1.350-0* 

13 

7091 

>9 

7094 

101*1371.223 

10141371.602 

0.379 

0.03k 

7.210-09 

19 

7091 

»*4 

7099 

10199643.124 

10199642.536 

-0.3*7 

0.023 

3 .**0-09 

40 

7091 

■4 

7119 

3929726.600. 

3929728.019 

-0.7*2 

0.03V 

2. 3*0-0* 

41 

7091 

• 4 

7119 

3 900996. 665 

3900597.570 


0.03* 

2.070-0* 

44 

7091 

*■» 

7120 

7660273.826 

7540273.123 

-0.701 

0.02V 

9.100-09 

43 

7091 

*•> 

79o7 

6267037.782 

6287020.271 

-17.311 

u.Ou* 

3.330-0* 

44 

7091 

*4 

79-t3 

12269596.212 

12249596.272 

0.05V 

0.02* 

*. 3*0-09 

43 

7094 

•4 

7099 

3616666.68b 

3514554.371 

-2.3ko 

0.02/ 

U**0-0* 

24 

7092 

»4 

7il«t 

7679017.59b 

7479018.461 

0.41*3 

0.02? 

*.300-09 

21 

7094 

*•0 

7113 

7686680.410 

7584661.155 

0.7*5 

0.02* 

*.7*0-09 

49 

7094 


7l*.U 

4015538.430 

4015538.979 

0.330 

0.0214 

1.3*0-08 

29 

7092 

»4 

790 f 

11 17 1115.715 

11171110.424 

-3.29k 

0.0*1 

*.7*0-09 

30 

7o94 

• 4 

7993 

5192643.026 

5192640.982 

-2.0V» 

0»o2* 

1 * 1*0-0* 

31 

7094 

*4 

7114 

74146Vb.9Sl 

7414696.912 

-0.0*0 

0.023 

7. **0-09 

34 

7094 

»*> 

7113 

7402692.901 

7402692.731 

-0.170 

0.0*3 

*•030-09 

33 

7094 

• 4 

7140 

4112220.642 

4112220.461 

-0.0*1 

0.023 

1 .**0-08 

34 

7094 

**4 

7907 

4373096.052 

9373093.497 

-0.35* 

0.0** 

1 .110-0* 

33 

7094 

«4 

7993 

4 554571.701 

4554572.165 

0.*** 

0.02* 

1, 2*0-0* 

34 

7114 

<m4 

7113 

256289.956 

266290.167 

0.210 

0.03* 

3. 330-07 

37 

7114 

'«k«4 

7140 

4022969.627 

4022959.505 

-0.022 

0.03k 

l.*ll>-0* 

34 

7114 

i4 

7907 

7 2*3 602.178 

7263566.024 

— I*.k3* 

0.U76 

2. *90-0* 

39 

7 114 

*4 

7993 

10667702.281 

1056 1702.70.1 

U ■ *20 

O.uU 

*.030-09 

40 

7113 

«*> 

71*0 

4096906.174 

6096V06.U6 

-0.0*7 

0.031 

1.790-0* 

4k , 

7U3 

■9 

7907 

703*726.657 

7036712.246 

-i*.*U 

O.of* 

2.310-0* 

42 

7113 

•9 

7943 

10596990.172 

10595990.425 

0.233 

0.01* 

h .070-09 

43 

7140 


79ul 

9047607. but 

9097399.369 

-*.2i2 

0.05* 

1.460-0* 

44 

7140 

>9 

7993 

7660V88.6V9 

7660969.3U0 

O.tOi 

0.022 

0.730-uV 

43 

7907 

>4 

7993 

19787493.056 

10767496.735 

4.616 

0.0*1 

*.9*0-09 
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Table 10 Residual Summaries by Station 


( a ) CONSOLIDATED STATISTICS FOR STATION * 706] 

PASS OBSEHV RESID MEAN RMS DEVIATION LENGTH MIN RfiSO MAX RFSO MEAN CLOS 


1 

4 

*0.9200 

5.30 i 

6.123 

1292.00 

*-8.493 

5.007 

-0.93 

2 

477 

0.1070 

0./137 

0.212 

1641.00 

-2.866 

0.398 

0.11 

1 

202 

0.1237 

0. 64 9 

0.632 

1494*00 

-6.840 

5.480 

0.12 

4 

6 

*2.0976 

3.459 

3.012 

1689.00 

-5.903 

1.400 

-2,10 

5 

859 

0.1242 

0.225 

0.187 

2358.00 

-2.436 

0.473 

0.12 

6 

1 

0.0450 

0.046 

0.0 

0.0 

0.046 

0.046 

sn 

o 

. 

© 

7 

1550 

0.0119 * 

0.322 

0.321 

2B10.00 

-4.383 

8.907 

0.01 

0 

4 

-4.4022 

5.625 

4.043 

1503.00 

-9.652 

-1.045 

i 

*• 

. 

© 

9 

14 

-0.4902 

2.473 

2.514 

2550.00 

-4.545 

5.946 

-0.50 

10 

1167 

-0.1706 

0.464 

0.4J2 

2484.00 

-6.694 

7.124 

-0.17 


(b) 

CONSOLIDATED STATISTICS 

FOK STATION 

: 7090 





PASS 

OBSERV 

RESID MEAN RMS 

DEVIATION 

waa.a.u.aus 

LENGTH 

MIN RESD MAX RCSO 

MEAN CLOS 

l 

97 

0.0802 

0.130 

0.095 

2156.00 

-0.221 

0.344 

0.09 

2 

167 

-0.0325 

0.104 

0.099 

2182.00 

-0.337 

0.177 

-0.03 

3 

102 

-0.0092 

0.131 

0.096 

2519.00 

-0.513 

0.131 

1 

o 

. 

o 

*0 

4 

207 

0.0322 

0.109 

0.105 

2471.00 

-0.282 

0.264 

0.03 

5 

196 

0.0335 

0.140 

0.137 

2631.00 

-0.433 

0.440 

0.03 

6 

141 

-0.0032 

0,142 

0.116 

1810.00 

-0.386 

0.192 

-0.08 

7 

263 

-0.Q764 

0.119 

0.091 

2853.00 

-0.427 

0.177 

-0.00 

0 

119 

0.0891 

0.138 

0. 106 

2037.01 

-0.153 

0.304 

0«09 

9 

67 

-0.0531 

0.093 

0.070 

1040.00 

-0.235 

0.145 

i 

o 

• 

o 

VI 

10 

171 

-0.0780 

0.122 

0.094 

2787.00 

-0.471 

0.101 

-o.oa 

11 

136 

-0.0981 

0.138 

0.097 

2036.00 

-0.451 

0.101 

-0.10 

12 

203 

-0.0640 

0.112 

0.092 

2481.00 

-0.526 

0.145 

-0.06 

13 

50 

0.0123 

0.000 

0.079 

574, 00 

-0.157 

0.216 

0.01 

14 

29 

0.0674 

0.143 

0.128 

1287.00 

-0.219 

0.286 

0.07 

15 

136 

-0.0940 

0.124 

0.081 

2339.00 

-0.347 

0.092 

-0.09 

16 

104 

0.1106 

0.435 

0.422 

2645.00 

-3.875 

0.378 

O.ll 

17 

55 

0.1305 

0.233 

0.194 

1318.01 

-0.282 

0.690 

0.13 

la 

’ 162 

-0.1606 

0.204 

0.126 

227?. 00 

-0.472 

0.071 

-0.16 

19 

173 

0.0735 

0.136 

0.114 

2491.00 

-0.272 

0.326 

0.07 

20 

136 

-0.1598 

0.212 

0.140 

2409.00 

-0.567 

O.lll 

-0.16 

21 

155 

0.0943 

0.185 

0.159 

2421.00 

-0.323 

0.378 

0.09 

22 

41 

0.0313 

0.069 

0.062 

1090.00 

-0.113 

0. 143 

0.03 

23 

88 

0.1992 

0.218 

0.088 

1438.00 

-0. 102 

0.348 

0.20 

24 

233 

-0.0433 

0.112 

0.104 

2842.00 

-0.597 

0.195 

-0.04 

25 

115 

0.1264 

0.1,56 

0.092 

1947.00 

-0.168 

0.345 

0.13 

26 

6 

0.1585 

0. 172 

0.073 

93.00 

0.071 

0.230 

0.16 

27 

* 189 

0.0197 

0.120 

0.126 

2630.00 

-0.289 

0.294 

0.02 

28 

154 

0,008 7 

0. 123 

0.085 

2109.00 

-0.164 

0. 328 

0.09 

29 

154 

-0.0249 

0.223 

0.223 

2088.00 

-0.567 

0.398 

-0.02 

30 

214 

0.0590 

0.120 

0.114 

2700.00 

-0.379 

0.292 

0.06 


42 
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Table 10 (cont'd) 


(c) CONSOL IDATEO STATISTICS FOR STATION : 7091 

.MSS OBSERV RESIO MEAN RMS DEVIATION LENGTH MIN RESO MAX RESO MEAN CLOS 


t * 

137 

0.0824 

0.271 

0.259 

1145.00 

-0.750 

0.683 

0.08 

2 

352 

-0.0646 

0.182 

0.171 

1718.00 

-1.040 

0.346 

-0.06 

3 

240 

0.0505 

0.169 

0.162 

1879.00 

-0.545 

0. 450 

0.05 

4 

439 

-0.0202 

0.314 

0.313 

1965.01 

-1.142 

4.742 

-0.02 

(d) 










CONSOLIDATED STATISTICS 

FOR STATION 

t 7092 





PASS 

OBSERV 

RESIO MEAN RMS DEVIATION 

LENGTH 

MIN RESO MAX KESD MFAN CLOS 

1 

322 

0.2185 

0.263 

0.156 

1761.01 

-0.258 

0.527 

0.22 

2 

286 

-0.1712 

0.268 

0.206 

1153.00 

-0.845 

0.666 

-0.17 

3 

1273 

-0.0004 

0.239 

0.239 

2380.00 

-1.234 

0.986 

i 

o 

. 

o 

© 

4 

363 

-0.0324 

0.304 

0.303 

2271.00 

-1.293 

0.734 

-0.03 

5 

9 

0.0926 

0.331 

0.337 

174.00 

-0.367 

0*5*6 

0.09 


(e) 

CONSOLIDATED 

STATISTICS 

FOR STATION 

: 7096 





PASS 

UBSERV RESIO 

MEAN RMS DEVIATION 

LENGTH MIN RESl) MAX RESO MEAN CLCS 

1 

969 

0.0078 

0.169 

0.189 

2389. 00 

-0.583 

0.546 

0.01 

2 

461 

0.0359 

0.150 

0.146 

2008.99 

-0.731 

0.331 

0.04 

3 

268 

-0.1355 

0.257 

0.219 

1109.01 

-0.931 

0.313 

-0.14 

4 

91 

-0.3075 

0.391 

0.244 

652.00 

-0.953 

0.122 

. 

0 

1 

5 

45 

0.0547 

0.166 

0.158 

924.00 

-O.A51 

0.356 

0.06 

6 

616 

0.0351 

0.213 

0.210 

1368.01 

-1.019 

0.531 

0.04 


(f) 

CONSOL IDATEO STATISTICS 

FUR STATION 

s 7114 





r^PASS ■ 

OBSERV 

RES ID MEAN 

RMS 

DEVIATION 

-,-iSNCTH 

MIN RESO 

t MAX t RE_SD 

MEAN CLOS 

1 

182 

-0.0405 

0. 176 

0.172 

1181.99 

-0.408 

1*00.3 

-0.04 

2 

17 

-0.1346 

1.475 

1.514 

1136.00 

-4.979 

2.632 

«n 

—A 

. 

0 

1 

3 

855 

0.0155 

0.263 

0.262 

2535.00 

-3.965 

1.966 

0.02 

4 

9 

1.1201 

2.490 

2.358 

1367.00 

-1.729 

5.392 

1.12 

5 

161 

0.0939 

0.155 

0.124 

1009.99 

-0.310 

0.465 

0.09 

6 , 

390 

-0.0036 

0.129 

0.129 

2102.99 

-0.3B7 

0.968 

o 

o 

• 

0 

1 

7 

6 

-0.7358 

3.655 

3.922 

88 7.00 

-4.838 

6.072 

-0.74 

8 

228 

-0*0292 

0.350 

0.349 

1045.00 

-0.384 

4.296 

-0.03 

9 

7 

0.0771 

0.111 

0.086 

676.00 

0.00 L 

0.236 

o.oa 

10 

7 

0.4243 

0.703 

0.606 

1080.00 

0.099 

1.779 

0.42 

11 

4 

-4.5052 

5.795 

4.208 

400.00 

-8.096 

0.651 

-4.51 
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Table 10 (cont'd) 


(g) 

CONSOLIDATED STATISTICS 

FOR STATION 

» 7115 





MSS OBSERV 

RES ID MEAN 

» a aaaaaa 

R^S DEVIATION LENGTH MIN RESO MAX RfSD MEAN CIOS 

1 

26* 

0*0850 

0*1)1 

0.099 

122*. 01 

-0.1/8 

0.379 

0.08 

2 

29 

0*2660 

0*280 

0*088 

711.00 

0.077 

o.*io 

0.27 

J 

364 

-0*0730 

1.057 

1.056 

2001.00 

-6.960 

0.595 

-0.07 

4 

27 

-6.7929 

6.793 

0.090 

1021.00 

-6.969 

-6.617 

-6.79 

— 

500 

0*0934 

0*151 

0.119 

2271.00 

-0.422 

1.167 

0.09 

* 

38 

0.1169 

0*328 

0.088 

1313.00 

0.120 

0.468 

0.32 

7 

172 

0.3589 

0*511 

0.364 

13*6.00 

-0.525 

*.82* 

0.36 

8 

63 

-0.0*68 

0*145 

0*138 

889.00 

-0.346 

0.257 

i 

© 

. 

© 

vn 

9 

119 

-0.2538 

0*312 

0*181 

1913.00 

-0.739 

0.504 

-0.25 

10 

588 

0.1T08 

0*213 

0.128 

2727.00 

-0.255 

1.611 

0.17 

11 

37 

0.1J03 

0.165 

0.102 

1608.00 

-0.075 

0.363 

0.13 

12 

44 

0.2*36 

0.286 

0*151 

652.00 

-0.040 

0.949 

0.24 


(h) 

CONSOLIDATED 

STATISTICS 

FOR STATION 

s 7120 





PASS 

M.HM 

QOSERV RES 1 0 MEAN 

RMS 

DEVIATION 

LENGTH MIN RtSD MAX RfcSD MEAN CLOS 

iimMaa«iaasj»aa.miMiiaiiMaiisii.iii 

1 

225 

-0.1213 

0* 142 

0.073 

1964.00 

-0.380 

0.120 

-0.12 

2 

44 

-0.0037 

0*157 

0.159 

829.00 

-0.772 

0.180 

-0.00 

3 

160 

0.0996 

0*140 

0*098 

1297.00 

-0*225 

0.313 

0.10 

4 

42 

-0*0689 

0*106 

0.084 

616.QQ 

-0-218 

0.698 

-0.07 

5 

IB 7 

0.0268 

0*133 

0* 13 1 

2614.00 

-0.3J7 

0.657 

0.03 

6 

346 

-0.0766 

0*114 

0.085 

2759.00 

-0.348 

0.247 

-0.06 

7 

401 

0.093 1 

0*138 

0.102 

2573.00 

-0.259 

0*298 

0*09 

a 

50 

-0.2583 

0.294 

0*141 

865.00 

-0.511 

0.034 

-0.26 

9 

121 

-0.1879 

0.216 

0*107 

1655.00 

-0.380 

0*104 

-* 

« 

0 

1 

10 

328 

0*1102 

0*163 

0*120 

2417.00 

-0.221 

0.409 

0.11 


(i) 

CONSOLIDATED STATISTICS 

FOR STATION 

: 7907 





PASS 

OBSERV 

RESIO MEAN RMS 

DEVIATION 

LENGTH 

MIN RESD 

MAX RESO 

MEAN CLOS 

1 

12 

0.1231 

0.366 

0*360 

292.40 

-0.406 

0*939 

0.12 

2 

41 

-0.0135 

0*604 

0*611 

1132*95 

-1.520 

0.761 

-0*01 

3 

51 

0.09L 7 

0*509 

0.505 

662.51 

-1.291 

0.603 

0.09 

4 

19 

-0*0486 

0.335 

0.341 

1027.76 

-0.767 

0.561 

-0.05 

5 

5 

-0*0084 

0*596 

0.666 

360.15 

-1.164 

0*446 

-0.01 

6 

52 

0*1667 

0*521 

0.498 

892.90 

-0.951 

2.032 

0.17 

7 

24 

0*0601 

0.401 

0.405 

660.27 

-0.925 

0.724 

0.06 

8 

19 

-0*2028 

0*393 

0.346 

607. 64 

-0.868 

0.360 

-0.20 

9 

35 

0*0493 

0*288 

0.287 

1162.58 

-0.648 

0.640 

0.05 

10 

53 

-0*0616 

0.408 

0*408 

914.95 

-1.008 

0.881 

-0.06 

U 

34 

-0.0226 

0*398 

0.403 

967.53 

-1.199 

0.697 

-0.02 

12 

5 

0*2501 

0*257 

0.065 

360.10 

0.169 

0.303 

0*25 

13 

24 

-0.5918 

1*005 

0.830 

1110.02 

-2.043 

0.962 

■m 

O 

1 

14 

17 

0*0004 

0*813 

0*638 

689.98 

-1.536 

1.569 

0.00 

15 

8 

0,2717 

0*526 

0.484 

420.00 

-0.434 

1.126 

0.27 

16 

9 

-0.U24 

0*568 

0*590 

465.05 

-1.312 

0.452 

-0.11 

17 

22 

-0.1859 

0*744 

0.737 

847.96 

-1.970 

1.663 

-0.19 

18 

28 

0.2387 

0.405 

0*333 

1012.50 

-0.646 

0*796 

0.2 V 

19 

2 

-0.3987 

0*425 

0*207 

90.00 

-0*545 

-0*252 

o 

• 

0 

1 

20 

29 

0*0272 

0*411 

0.417 

779.97 

-0.855 

0.698 

0.03 
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Table 10 (cont'd) 
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CONSOL I OAT ED 

</* 
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•4 

O 

(A 

FOR STATION 

x 7943 





PASS 

08SERV RES 1 0 MEAN 

RMS 

DEVIATION 

LENGTH 

MIN RESO 

MAX RESO 

ME At* CLUS 

mmmmmmm a 









l 

56 

0.0012 

0.433 

0.437 

1154.92 

-1.019 

0.959 

0.00 

2 

122 

0.0426 

0.338 

0.336 

2325. 29 

-1.432 

0.815 

0.04 

3 

82 

-0.0668 

0.439 

0.437 

1905.10 

-1.252 

t.074 

-0.07 

.4 

98 

-0.0301 

0.408 

0.409 

2272.34 

-0.970 

0.715 

-0.03 

3 

83 

0.0633 

0.426 

0.424 

1484.75 

-1.229 

1.222 

0.06 

4 

85 

0.0516 

0.337 

0.335 

1547.79 

-0.917 

0.872 

0.05 

7 

90 

0.1375 

0.292 

0.259 

1424.79 

-0.489 

0.99! 

0.14 

• 

212 

-0.0796 

0.318 

0.309 

2744.95 

-0.564 

0.679 

-0.08 

9 

159 

0.0983 

0.344 

0.330 

2910.00 

-0.720 

0.861 

0.10 

10 

47 

-0.0243 

0.511 

0.516 

1934.86 

-1.052 

1.198 

-0.02 

a 

54 

0.0337 

0.429 

0.432 

1514.65 

-1.193 

1.353 

0.03 

12 

75 

-0.1024 

0.552 

0.546 

1605.00 

-1.306 

0.974 

-0. 10 

13 

139 

0.0414 

0.239 

0.236 

1822.48 

-0.681 

0.594 

0.04 

14 

100 

-0.0034 

0.375 

0.377 

1897.49 

-0.851 

1.041 

-0.00 

IS 

170 

-0.0420 

0.414 

0.413 

2955.00 

-1.097 

1.295 

-0.04 

16 

106 

-0.0837 

0.347 

0.338 

1890.00 

-1.137 

0.715 

-0.08 

17 

90 

-0*0132 

0.280 

0.291 

1919.92 

-0*878 

0. 702 

-0.01 

i$ 

44 

-0.0889 

0.470 

0.467 

2167.47 

-1.129 

0.643 

-0.09 

19 

79 

0.1667 

0.424 

0.392 

2287.55 

-0.943 

0.899 

0.17 

20 

178 

0.02L2 

0.282 

0.282 

2489.80 

-0.607 

1.253 

0.02 

21 

30 

-0.1376 

0.625 

0.617 

1492.55 

-1.638 

1.460 

-0.14 

22 

88 

-0.0665 

0.485 

0.483 

2092.54 

-1. 147 

1.249 

-0.07 

23 

172 

-0.0355 

0.263 

0.261 

2827.32 

-0.o7l 

0.611 

-0.04 

24 

73 

-0.3053 

0.399 

0.259 

629.94 

-0.096 

0.119 

-0.31 

25 

67 

0.0351 

0.448 

0.450 

1619.87 

-1.124 

0.798 

0.04 

26 

96 

-0.0945 

0.596 

0.592 

1619.69 

-1.605 

1.349 

-0.09 

27 

163 

-0.1419 

0.279 

0.241 

1919.90 

-0.771 

0.511 

-0.14 

2$ 

191 

0.1105 

0.330 

0. J 12 

1994.97 

-0.72 7 

0.898 

0.11 

29 

262 

0.1799 

0.365 

0.U8 

2999.78 

-0.819 

0.837 

0.18 

30 

97 

0.1417 

0.282 

0.245 

1942.65 

-0.673 

0.703 

0.14 

31 

61 

0.1398 

0 • 466 

0.448 

2175.10 

-0.96 1 

0.875 

0.14 

32 

41 

-0.0845 

0.405 

0.483 

1507.4V 

-1.351 

0.607 

-0.09 
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2.22 Doppler Experiments 


A 3-*/3 

' 1*83 13540 


2.221 Geometric Adjustment of Simultaneous Doppler-Derived 
Range Differences 


The results of work on this topic are described in a paper presented 
at the Third Internationl Symposium on the Use of Artificial Satellites 
for Geodesy and Geodynamics, Ermioni, Greece, September 20-25, 1982. It 
appears on the following pages and will be published in the proceedings 
of the symposium obtainable from the National Technical University, Athens. 


46 


TkOid TnteAncutionaZ Symposium on the. Uie o<{ A>Ui6icMit SatoZLitzA ion 
GtodzMj and GcodyncurUcb , Efmioni, Greece, SzptmozA 20-25, 1962 


GEOMETRIC ADJUSTMENT OF SIMULTANEOUS DOPPLER-DERIVED 

RANGE DIFFERENCES 


^Chengze Zhang* and Ivan I. Mueller 
Department of Geodetic Science and Surveying 
The Ohio State University, Columbus, Ohio 43210 USA 


ABSTRACT. A mathematical model for the use of simultaneous Doppler- 
derived correlated ranges in the geometric mode is presented. The model is 
tested with data taken during the EDOC-2 campaign with different integration 
intervals. The results of this adjustment are compared with the EDOC-2 
adopted solution and those from an uncorrelated model [Schneeberger et al< 
1982] used earlier to provide more economical calculations. 

The analysis of the comparison shows that the correlated mode is 
superior to the uncorrelated one when the optimum integration interval of 
23 seconds is used. 


1, INTRODUCTION 

The geometric purpose of satellite geodesy is to tie remote stations 
together in the same geometric system. Its ultimate aim is to determine 
the coordinates of unknown ground stations [Mueller 1984]. 

Satellite geodesy with Doppler techniques is based on the principle 
that a frequency transmitted from a satellite-borne transmitter moving 
relative to a ground receiver is observed shifted by the Doppler effect. 
The observations are Doppler counts which are measures of the range change 
between the satellite and the receiver during the integration interval 
[Wells 1974]. 

In the geometric mode for Doppler observations, the satellite is 
regarded as a benchmark in space and its coordinates at the observation 
instants are unknowns which are solved in an adjustment with the unknown 
coordinates of ground stations. Such solutions are based on geometric 
rather than dynamic principles; therefore the calculations are relatively 
simple and do not require extensive computer programs. 


*0n leave from the Institute of Geodesy and Geophysics, Chinese Academy of 
Sciences. Wuhan. People's Republic of China. 


In a previous study [Schneeberger et al. 1982], the Doppler-derived 
ranges were regarded as uncorrelated pseudo-observations as a further 
simplification (to save computer time). In fact, since the Doppler-derived 
ranges are calculated from Doppler counts, it is obvious that there exist 
correlations between them in a given pass. The purpose of this study is 
to investigate the use of Doppler-derived correlated ranges in the geometric 
mode. 


This method is then tested against a data set from which a dynamic 
solution is available. The results are compared with both the dynamic 
solutions and the uncorrelated geometric one. 


2. SUMMARY OF THE PREVIOUS STUDY BASED ON UNCORRELATED OBSERVATIONS 
[SCHNEEBERGER ET AL. 1982] 

2.1 Definitions 

The coordinate, tyttm in which the computations are performed is an 
earth-fixed Cartesian system. It is defined by the assigned six coordinates 
distributed among at least three ground stations. A iatMite point is the 
position of a satellite at a certain epoch. An even* is the set of all obser- 
vations to a satellite point* A pea* is a set of satellite points between two 1 
epochs which are observed without Interruption from at least six ground sta- 
tions. A VoppieA- derived range is a pseudo-observation derived by adding the 
range differences computed from Doppler counts to an estimated initial range. 

2.2 Doppler-Derived Ranges 


The basic equation which related the ratio between the received frequen- 
cy f and the transmitted frequency fo to the range rate between transmitter 
and receiver (f) is accredited to Doppler (1803-1853); 


f 

77 


(rh) * (* - I) 


where c is the velocity of propagation for electromagnetic waves in a vacuum. 
This equation has to be integrated to find a relation between the shifted fre- 
quency and the range difference during a time interval t. A detailed derivation 
can be found in [Brown and Trotter 1969] resulting in 


r j - rj_j * X 0 (Nj - afoo tj) + S , (1) 

where 

r. = range between receiver and transmitter at epoch T. 

3 J 

r j-l “ range at e P 0C h 

N. B the Integrated Doppler shift over time interval tj = T* - T* , 

* (referred to as the Doppler count) J J 
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Afoo * the difference between the transmitted frequency and the reference 
frequency generated in the receiver 

■ I 4 * * wavelength corresponding to the frequency of transmission f® 
c 

S * correction term representing all systematic errors such as bias in 
the difference between the adopted transmitted and reference fre- 
quencies, and/or the drift rates of transmitter and receiver frequen- 
cies, 

Substituting the range difference computed from the Doppler count 
Ar^ * X 0 (Nj + Afoo tj) 


Into equation (1), the range at epoch Tj is 
r j • r M + ir j + s j 


If the range ro at an initial epoch To Is known, the range for an epoch T k can 
then be calculated from (taking into account that most Instruments reset the 
Doppler count for each Interval) 


r 0 


k 

% 


Ar * Si. 


( 2 ) 


This equation is correct only in a vaGuum. Since the signal is passing through 
the ionosphere and the troposphere, the range has to be corrected for refrac- * 
tive effects. The ionospheric refraction is automatically compensated (to 
first order) by measuring the Doppler shift of the two different frequencies 
(400 and 150 MHz) [Krakiwsky and Wells 1971]. Each range has to be corrected 
therefore only for the tropospheric refraction AT... The tropospheric refrac- 
tion model used in this study is the one outlined in [Brown and Trotter, 1973], 
using the Smith-Weintraub model for the index of refraction [Jordan et al. 
1966]. 


Since the initial range in equation (2) is not known, we must use- an 
approximate Initial range ro and add a correction term ao to be estimated from 
the adjustment, 


r 0 = ro + a 0 

a 0 is considered part of the systematic error term S k in equation (2). The 
modelling of the other systematic effects in Sk is , given in great detail in 
[Brown and Trotter 1969; Kouba and Boal 1976]. In this study only two major 
terms are used: a + bt. The main cause of the constant term a is the possible 

bias in the adopted frequency fo, and the initial range error a 0 above. The 
time dependent term bt is caused mainly by the difference in the adopted val- 
ues for the transmitter and receiver frequencies (frequency offset) Af 0 ®, 
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Other terms in the systematic error model mentioned by Brown and Trotter [1969] 
but not considered in this study are range dependency, a function of the 
second power of time, and a function of the elevation angle (for residual 
refraction errors). An explanation of why only the above two terms are used 
here may be found in [Schneeberger 1982]. 

Substituting all terms for S and the correction for tropospheric refrac- 
tion equation (2) can be written as 


Ik 


k 

ro + £ Ai'j + ATp + a,j + bjt^ 

j*l 


where the subscript 1 refers to ground station 1 . Defining the VopplttArdeALvtd 
Aange as 


k 

r 0k * r ® + I Ar j + AT r » (2‘) 

and recalling that 

and changing the signs of a and b, we arrive at the mathematical model 

«-D 1k . A\ - x,)‘ + (Y k - Y^)* + (z k - Z^) s f a u + b u t k (3) 

where r0j k Is the Doppler-derived pseudo-range (derived from the measured 
Doppler 1K counts and corrected for tropospheric refraction), and the unknown 
parameters to be solved for in a least squares adjustment are 

Xj, Yj, Z 1 the’ unknown station (1) coordinates 

fcjj, Y k , Z k the unknown satellite (k) coordinates 

a.,, b,. . the unknown coefficients used to model systematic errors 

u for each station (1) and pass (a) 

t k is the time elapsed from the epoch of the initial range r<». 


2.3 Least Squares Adjustment 

The mathematical model developed above has the form of an observation 
equation: 

L a = F(X a ) (4) 


ORIGINAL PAGE IS 
0 ¥ POOR QUAunr 
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where L Is the adjusted Doppler-derived range, and X is the vector of the 
unknown parameters which can be divided into three a subvectors: 


XG a » XG 0 + XG 
xc a » XC 0 + XC 
XS a - XS 0 + XS 

a 


containing the coordinates of the ground stations 
containing the error coefficients a, b 
containing the satellite coordinates 


Equation (3) can be written in linearized form 


r Dik* + v ikt ■ F ikt + Ik v* xg i + m v' xc u + Ire * xs k * ••• 


*Xo 


Xo 


'Xo 


(5) 


or, after neglecting higher-order terms 

v jka “ A ikJt' XG 1 + C 1 kJt* XC i A + S ikA* XS k " W ikA 


where 


R iki 3XG 


xg? xc} A xsjj 


AM. -IfcIL -M\ 

' r "iU r Mkt r, 1kt ' 


c m a lfe| X6 . xc.^ xs . • (1 - V 


<. = 3F 
b iki are 


xg$ XC| A XS° 


-A 


1k£ 


M 1U * r Dikt " 1 ki + a °U + b °1]tV 

<\, jk = /fx“-xp< + (yj-vp 2 + (z»-z ») 2 


In this study all pseudo-range observations are assumed to have equal 
weight. For reason of convenience in programming, the a priori variance of 
unit weight is chosen to be equal to the variance of a range observation 

_ _2 

00 °DR 


Therefore all observations have the weight one. Further details of this 
least squares adjustment as used in the Geometric Doppler (GEODOR) computer 
program may be found in [Schneeberger 1982], 
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3. ADJUSTMENT WITH CORRELATIONS CONSIDERED 
3.1 Mathematical Model 


ORIGINAL PAGE IS 
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The correlation existing in the Doppler-derived ranges are considered 
in this study by assuming that the range differences (computed from the 
Doppler counts) are independent observations. 

Under this consideration, substituting eq. (2‘) into (3) and moving 
all the terms to the left side, we obtain 

/ <V X 1> 2 + <YV 2 + <VV 2 + «H + b u‘k ’ r " ‘ } ml 4r J - aT r - 0 

( 6 ) 

Thus the model becomes the form of a condition equation with parameters: 
F(l a . X„) = 0 

Eq. (6) can be written in a linearized form, using the same notation as 
before, 

A 1 U X 1 * C 1 k XC 1 k* + S ik<t XS k + j. B ik« - W ik* * 0 < 8 > 

J 1 


(7) 


where B..^ stands for the derivatives of F with respect to Ar 4 , i.e., 


B ila 


3F 


3Ar . 

J 


_ M if j < 

\ 0 if j > 


<. k 
k 


(9) 


All the observations are assumed to have equal weight. For convenience 
in programming, the a priori variance of unit weight is chosen to be equal to 
the variance of a range difference observation 

_2 _ _2 
so - o Ar 

Therefore, all the observations have unit weight. For the detail of the 
derivations of the mathematical model and the method of solving this prob- 
lem, see [Zhang 1982]. 

3.2 Construction of Normal Equations 

The solution of the normal equation system for the least squares 
model of condition equations with parameters has the following form 
[Uotila 1976]: 


where 


X = -(A T m" 1 A)"^ a t m -1 w 


M" 1 = (B P" 1 B 1 )" 1 


( 10 ) 
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Therefore, before constructing the normal equation system, M" 1 has to be 
found first. Fortunately, the matrix B has a regular configuration, and 
so does the matrix M" 1 [Ashkenazi et al . 1980]. For the sake of simplicity, 
we investigate a matrix B for one station and one pass. From eqs. (8) and 
(9) it is evident that the matrix B has the form 


B 


-1 0 

-1 -1 

-1 -1 

• * 

* ♦ 

, • 

-1 -1 


0 0 

0 0 

-1 0 

♦ * 

• i 

• • 

-1 -1 


... 0 

. . . 0 

. . . 0 

« • 

• • 

♦ 

... -1 


0 - 2 ) 


If one assumes uniform weight and no correlations between the range differ- 
ences, and chooses the variance of unit weight equal to the variance of 
range difference observation, the matrix P will become an identity matrix. 
Then the matrix M can be written as 


M = B P-V 


1 1 1 1 ... 1 

1 2 2 2 ... 2 

1 2 3 3 ... 3 



1 2 3 4 ... n 


(13) 


where n is the number of observations in this pass. is found by 
inverting M: 

I 2 -1 0 0 ... 0 0 0 I 



04 ) 


Since M is a regular diagonal matrix, it will not invite much difficulty 
when constructing normal equations. For the case of more than one station 
and more than one pass, matrices B and M * can easily be found by using the 
same method [Zhang 1982]. 


' -l 

After the matrix M is found, all the coefficients of the normal 
equation system can be calculated. Since this normal equation system is 
still of the sparsity pattern, a method called second-order partitioned 
regression can be used to eliminate the unknowns to save storage and 
computing time [Brown and Trotter 1969]. 
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4. NUMERICAL TEST 

4.1 Solutions and Their Comparisons 

The data taken during EDOC-2 was used for testing the uncorrelated 
and correlated modes. Fig. 1 shows the network used which is chosen from 
EDOC-2. There were many solutions for each mode, but only the best , 
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Fiq. 1 EDOC-2 network [Boucher et al . 1981] 


one of each mode can be presented here. Table lisa summary of these 
two solutions. Solution P4-5 is in the uncorrelated mode; solution C-5 
in the correlated mode. The integration Intervals of both solutions are 
5 X 4.6 * 23 seconds, 

The information of solutions using different integration Intervals 
from the correlated mode is collected in Table 2. In the designation 
C-1 , i indicates the Integration IfiM/vals used, e.g., in case of i * 2, 
the range change is over 2 x 4.6 = 9.2 s. Fig. 2 gives a visual comparison 
of these solutions. It is obvious that the solution with i * 5 is the 
best. 


4.2 Test of the Systematic Error Model 

In this study as in the earlier one only the two major terms are 
used for modeling the systematic effects: a + bt. The residuals of the 

observations of randomly selected passes from the total of 193 passes were 
plotted for each station. Fig. 3 is one example. Investigating the distri- 
butions of the residuals of the observations at each station, no signifi- 
cant remaining systematic effect is found, which indicates that the two 
major terms used for modeling the systematic effects are reasonable. 

4.3 Test of the Residuals 

From Table 1 we can find tnat the correlated mode is superior to the 
uncorrelated one. In spite of that, there are still significant differ- 
ences between solutions C-5 and EDOC-2. In order to find the reason, the 
residuals of all observations were investigated. Table 3 lists the 
statistics of the residuals of the observations over the ten worst passes. 

Checking this table, one can see that the maximum residual is as 
large as 160 m, and the ratio of the number of the observations whose 
absolute residuals are larger than three times the standard deviation, 
to the total number of the observations for each one of the worst passes 
is high. The worst one is as high as 11.2ft. This indicates that there 
may be blunders in the data set. 

4.4 Problem of Weights 

As stated earlier, all observations are assumed to have equal weight 
and the a priori variance of unit weight is chosen to be equal to the 
variance of a range difference observation 

Oo = <4 = 1.0 

In Table 2 one can see that the a posteriori standard deviations of unit 
weight for all the solutions are much larger than the chosen a priori one. 
For instance the a posteriori standard deviation of unit weight of the 
best solution, C-5, is as large as 3.4. 
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Table 1 Summary of Solutions F4-5 and C-5 


Solution No.: 

F4-5 

C-5 

Total No. of Passes. Processed 

193 


193 


total No. of Events 

3,430 


3,430 


No. of Unknowns 





Station Coordinates 

30 


30 


Error Coefficients 

3,312 


3,312 


Satellite Coordinates 

10,290 


10,290 


Total No. of Unknowns 

13,632 


13,632 


Total No. of Observations 

27,531 


27,531 


Degrees of Freedom 

13,899 


13,899 


A Priori Weight Information: 





Range (or Range Difference) 

1 m 


1 m 


Error Coefficient o a 

50 m 


50 m 


Error Coefficient 

38 rn / 2 min i 

38 m / 2 min 

Fixed Station Coordinates ox» oy, oz 

1 mm 


1 mm 


Other Station Coordinates o%, oy, oi 

100 m 


100 m 


3 Satellite Events/Pass ox» °Y> 

10 m 


10 m 


A Posteriori Standard Deviation 

3,5 m 


1. d m 


of Unit Weight 





Coordinate Differences with Respect 

AX 

aH 

A$ AX 

AH 

to EDOC-2 Solution (all units in m) 





Station No. 220* 

0.0 0.0 

0.0 

0.0 0.0 

0.0 

(* indicates fixed station) 221 

3.5 16.2 

-4.4 

3.0 7.0 

-4.6 

222 

3.4 -4.3 

-1.2 

1.9 -3.0 

-1.4 

223 

1.5 9.8 

-7.1 

1.7 5.2 

-4.2 

224 

5.6 27.1 - 

•12.6 

5.3 13.2 

13.1 

225 

-4.8 -10.9 

6.7 

-0.2 18.4 

3.0 

226* 

0.0 0.0 

0.0 

0.0 0.0 

0,0 

230 

1.0 -4.5 

13.6 

4.6 3.2 

9.9 

231* 

0.0 0.0 

0.0 

0.0 0.0 

0.0 

232 

4.5 -3.1 

15.7 

2.1 -4.6 

20.5 

233 

0.8 8.7 

-2.0 

1.6 2.0 

-0.9 

234 

0.1 -1.7 

5.3 

-1.0 1.1 

7.5 

235 

1.0 44.2 

70.0 

8.2 -'0.4 

1.9 

Average absolute difference (m) 

2.6 13.0 

14.0 

2.2 5.8 

6.7 

(10 stations) 

±2.8 ±5.2 

±4.6 

±4.7 ±8.9 

±7.9 

Average absolute difference 
In position (m) 

20.8 ±21.7 

10.8 ±6.6 

Average absolute statlon-to-statlon 

10 2 ±10 5 

S.5 ±4 5 


chord distance difference (m) 

lUib iJlU . W 
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Table 2 Comparison of the Different Integration Intervals Used in 
Adjustment 


Name of Solution: 

C-2 

C-5 

C-10 

C-15 

Integration interval (seconds) 

9.2 

23 

46 

69 

Computing time (minutes)* 

25,20 

8.83 

5.96 

2.79 

A posteriori standard deviation 
of unit weight 

2.4 

3.4 

4.5 

5.9 

Average total absolute difference 
in position (m) (10 stations) 

10.9 

±7.3 

10.9 

±6.6 

22.0 

±19.3 

33.8 

±37.1 

Average absolute station-to-station 
chord distance difference (m) 

6.4 

±5.6 

5.5 

±4.4 

8.9 

±8.2 

17.1 

±19.2 


c : computing time used 

B: average of absolute difference 

of baselines 




46 69 

integration interval (s) 


Fig. 2 Computing time used and average of absolute differences of 
baselines plotted against the length of integration interval 
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Distribution of residuals at station 223 in passes 4, 50, 101 and 181. 



Table 3 Statistics of the Residuals of the Observations of the 
. Ten Worst Passes 


Pass Number of Observations 


NO. 

No. 

Total 

v > 

3o 


V 

> 2a 

v > 

10.0 m 

^max 

Number 

% 

Number % 

Number 

'% 

l 

49 

207 

20 

9.7 

27 

13.0 

17 

8.2 

160.8 

2 

46 

187 

18 

9.6 

38 

20.3 

13 

7.0 

126.1 

3 

187 

143 

16 

11.2 

27 

18.9 

12 

8.4 

41.4 

4 

21 

262 

19 

7.3 

32 

11.8 

10 

3.8 

53.0 

5 

43 

181 

9 

5.0 

17 

9.4 

8 

4.4 

39.5 

6 

180 

142 

9 

6.3 

15 

10.6 

8 

5.6 

25.6 

7 

51 

221 

10 

4.5 

15 

6.8 

7 

3.2 

17.2 

8 

26 

186 

10 

5.4 

12 

6.5 

5 

2.7 

26.0 

9 

25 

105 

6 

5.7 


7 

6.7 

5 

4.8 

31.6 

10 

16 

195 

7 

3.6 

10 

5.1 

4 

2.1 

37.4 


Table 4 presents the comparison of the weights of each station 
calculated from the residuals over all passes. The weights of the stations 
differ from each other for solutions C-5; the largest one is ninefold as 
large as the smallest one. When the ten worst passes are taken out, the 
weights are close to each other, and the a posteriori standard deviation 
of unit weight is decreased from 3.4 to 2.0. It is seen that the exist- 
ence of blunders is probably the most important detrimental factor in the 
solution. 

Unfortunately, neither taking out the ten worst passes nor repeating 
the computation with the different weights for each station improved the 
result. It is likely that although taking out the ten worst passes removed 
the major blunders, it also resulted in losing many useful observations. 


ssssss 
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Table 4 Comparisons of the Weights of Each Station 


Station 

All Passes Used 

W/o 10 Worst Passes 

No. 

No. of Obs. 

0 

mm 


a 

0 

P 

220 

1579 

2.08 

2.7 

1473 

1.73 

1.3 

221 

1668 

1.72 

4.0 

1606 

1.56 

1.6 

222 

2912 

3.30 

1.1 

2711 

1.33 

2.2 

223 

1862 

1.30 

6.9 

1777 

1.27 

2.4 

224 

2821 

1.77 

3.7 

2631 

1.49 

1.7 

225 

1893 

2.19 

2.4 

1757 

1.33 

2.2 

226 

845 

2.56 

1.8 

763 

1.26 

2.4 

230 

2505 

2.26 

2.3 

2435 

1.21 

2.7 

231 

2789 

3.85 

0.8 

2609 

1.87 

1.1 

232 

2391 

2.42 

2.0 

2205 

1.30 

2.3 

233 

2760 

1.67 

4.2 

2575 

1.30 

2.3 

234 

2641 

2.36 

2.1 

2444 

0.99 

4.0 

235 

865 

1.25 

7.5 

809 

1.01 

3.8 

Degree of 
Freedom 

13,899 

12,910 

A 

Oo 

3.4 

2.0 


5. CONCLUSIONS 

On the basis of the comparisons, the folloiwng conclusions can be 

drawn : 

(1) The geometric mode of solving the problem of simultaneous Doppler 
derived ranges without considering the correlation is a weak one. 

(2) The correlated geometric mode leads to better results. Comparing 
with the uncorrelated solution, the correlated mode reduced the average 
total absolute differences (with respect to EDOC-2) in position from 20.8 
±21.7 m to 10.9 ±6.6 m; and the average absolute station-to-station chord 
distance differences from 10.2 ±10.5 m to 5.5 ± 4.5 m. 

(3) The choice of the optimum integration interval is very important 
for the use of simultaneous Doopler-derived ranges in the geometric mode. 

The examples of this study demonstrate that the optimum integration inter- 
val is 23 s, which agrees with that suggested by [Ashkenazi et al. 1980]. 
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2.222 Doppler Interccmpariscn Experiment 


In the previous semi-annual report, preliminary results 
of the 1979 CSU comparison test of Doppler receivers are 
given. Since that time a final report on the comparison has 
been completed [Archinal, 1982] as a Master's thesis (and 
soon as a report of the Department of Geodetic Science and 
Surveying) . 

In this report, seme of the results presented in the pre- 
vious report ace revised, and seme additional final results 
are presented as well* For a more detailed discussion of 
the following, refer to [Archinal, 1982], and [Acchioai and 
Mueller, 1982]. 


AL_RESU LI S_CF_DAj[ A__ SEDUCTION 

As mentioned above, some of the results presented here 
are slightly different than those given in the last report. 
This is primarily due to: 

a) The determination and use of receiver time delays in 
the GEODOP processing. 

b) The modification of GEODGP to allow the input of a 
"common station noise" estimate, and use of this op- 
tion, along with the use of a variance estimation pro- 
cess in GEGDCE as well. 


Therefore revised versions of the tracking statistics and 
chord difference results are given here, along with new ia- 
foraatior concerning the estimation of the receivers' obser- 
vational (range rate) error and oscillator stability. 
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3 £3£JsAm_ statists 

The statistics on the number of passes tracked# used in 
P8ED0P# and two types of GEODOP runs are presented in tacle 
t. Although several numbers have changed substantially from 
those given in the last report# most of the results given 
there are still valid. In addition to these results# it 
should be noted that if these statistics are broken down by 
antenna setup (as in [Archinal# 1982# pp, t>1-70]), it be- 
comes clear that: 

a) The CMA-751 and the JMB-lAs generally tracked about 
the same number of passes# and slightly acre than the 
MX1502 (when operating correctly and tracking continu- 
ously). 

b) Ibere is no bias due to antenna location# at least 
when PiEHOP rejections are considered. Tne relative 
percentages cf rejections stayed fairly constant for 
all setups for the JMB-1A #2 and MX15G2. No conclu- 
sions can be drawn for the C HA- 751 due to its faulty 
antenna cable (cn all but one site)# or for the JMB-1A 
HI# since it only occupied oue site. 


The GEODCJ? statistics (for a multi-station broaucast 
epheiaeris solution and a single station precise ephemeris 
solution) show a fairly ccnsistant observation*’ miss value 
for all instruments# except for the MX 1502, wtiich has a 
higher value in both solutiens. This higher value is due to 
the fact that the MX1502 was recording only (the tetter) 
passes which Reached over 15 degrees altitude cn the first 
setup# which strongly affects the grand totals showu here. 
The observations/pass for the CMA-751 are not representative 
here either# since it was operating properly only during the 
first and last setup. 


Ch or d Di. ffere nc e_Besu j .ts 

Table 2 shows absolute differences obtained in the chord 
distance between all pairs of instruments for each antenna 
setup for multi-station and precise ephemeris solutiens. 
Many of these values are different from those given in the 
previous report, with generally smaller standard deviations 
and chord differences than previously reported. This is 
probably due to the changes in weighting and the tetter de - 
cermined delays respectively, and points out the value cf the 
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more riyorus weighting and tetter determination of delay for 
these solutions. 

Even with the different* s, the previous result still 
holds, that most of the differences (all execept two) lie 
within their three sigma value. A new result is that the 
single baseline determined between instruments of the same 
type (between the two JME-IAs on setup #1) did net show sig- 
nificantly better results than pairings between any ether 
instrument combination. In conclusion, it appears that 
there is no evidence that any of these instruments are bi- 
ased against one another for chord determinations (over 
short distances) . 

Another additional result shown by this table is that the 
precise epheiueris two satellite solutions for chord distanc- 
es do not appear to have necessarily higher accuracy cr pre- 
cision than the corresponding broadcast ephemeris five sat- 
ellite solutions, and in fact the precision cf the broadcast 
ephemeris solution is better in all cases. This simply in- 
dicates that the greater number of observations in the 
broadcast ephemeris solution improves the results mere than 
the corresponding increase in ephemeris accuracy of the pre- 
cise ephemeris solution. This would imply that if cnly 
chord distances were needed from Doppler observations , then 
generally broadcast ephemeris solutions would be preferable 
to precise ephemeris solutions, since the former usually 
have more data available. 


Eange Ba te Meas urement E r rors 

Using procedures described in detail in [Archinal, 1982, 
pp„ 70-79], estimates of the common station noise and each 
instrument's range rate standard deviation were mads for 
each setup and precise ephemeris satellite. The results are 
shown in table 3 and discussed here. 

First of all, the common station noise was estimated by 
processing only simultaneous observatins and precise ephem- 
eris orbits. The common station (cr "interstation" cr "sat- 
ellite" noise) estimates were made using the common station 
estimated variance-covariance matrix output by GECEGE tc ob- 
tain the values shown in column three of table 3. The re- 
sults vary with satellite and time during the entire test, 
with an amount between 3.4 and 7.5 cm/30 seccnds. The ove- 
rall average value (weighted mean of all observation pairs) 
is 4.9 cm/30 seconds. Since the range was not tcc great. 
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and rather than change the value tor each setup/satellite 
(perhaps based on too few observations) , the GJSCCOP default 
of 5-0 cm/30 seconds was then used for all subsequent pro- 
cessing. 


Secondly, the estimated receiver range rate standard de- 
viations were computed for each setup and satellite, using 
the 'weighted mean of the diagonal elements cf the estimated 
variance-covariance matrix of the residuals, 
for all three instrument types are shown 
columns of table 3. These values were ob! 
single station precise epheraeris solutions, 
servations were approximately (the reject 
tical testing cause some exceptions) simi 
some of the individual solutions did not 
be considered significant, using at least 
in each case (but only 350 for the JHE-1A # 
range rate standard deviations were found 
10.4, and 12.5 cm/3u seconds for the C 
JMR-#2, and the dX1502 respectively, over 
of the test. The variations shown during 
partially instrument related, Lut they 
mostly to the satellite noise gust discussed, 
precision of the three instruments continuously 
also stays approximately the same for all 
satellites, which also indicates that tl 
non-receiver related. It is also significant 
ation between instruments is usually less 
onds, showing that these i.nstrss'gents ace c 
milar, and that the variation in the ccamoi 
generally greater than this. The conciui 
be drawn that the variation of the measurement precision 
between these instruments is not si gn if icant. The even more 
important conclusion which can be drawn is that the range 
rate accuracy obtainable depends in many cases more on the 
time and the satellite chan it does on the receiver itself. 
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Lastly, to obtain the best possible estimates cf the fi- 
nal variance-covariance matrices in GECDGP , the GECDCP op- 
tion was used to allow an internal estimate of the range 
rate standard deviation and adjacent observation correlation 
for each pass to be made, with the previously estimated 
range rate standard deviation value (given in the last para- 
graph) used as an input acpoximate value. Although increas- 
ing the computational time by over 50* (all of the passes 
are processed twice), this method takes into account the 
variation of the satellite noise and possible variations in 
the receiver noise during the period under consideration. 
It is felt that this procedure, in conjunction with the 
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first two above would result in the ucst rigcrus processing 
of the data, to provide the best solutions. 


lJ&2USaSX_££i£t_.&££Ui£s 

The frequency drift of an instrument’s oscillator is an 
important quantity which can be determined to fairly high 
accuracy during data reduction. In general, the more stable 
an oscillator (over the period of a satellite pass) the bet- 
ter the timing and Doppler count measurements can be made. 
If a drift is occurring, and remains fairly constant in 
time, it can bo taken into account m the adjustment of the 
data (as in GEODGE) , although it must still be assumed tc be 
linear over a pass, and should not be very large in magni- 
tude. If the drift is erratic, either changing during a 
pass or over gust a few passes, the data will be very noisy 
due to these unmodeled changes in the oscillator. It is 
therefore important to check the frequency drift variations 
of these instruments. Ideally, one would like to check the 
short term drift which corresponds to the length of a satel- 
lite pass (over 100 seconds to about 15 minutes), but this 
is generally not possible unless an atomic standard is 
available for comparison. Instead, the long terra drift of 
these instruments can be checked for variations (which may 
provide an indication of the short term stability), or at 
least checked against the manufacturer's specifications. 

In the case of the data collected here, the frequency 
drift for each instrument for eacb setup and precise ephem- 
ecis satellite has been determined. The values have been 
obtained from the difference between the first and last 
(reasonable) frequency offsets computed fci each instrument 
during a setup. Ihe frequency offsets were determined from 
two satellite (one satellite at a time) precise ephemeris, 
single station solutions, and the antenna setup periods 
which ranged from about five to fifteen days in length. 
Note that to ortaiu the per day values given here, the as- 
sumption has been made that the frequency drift is constant 
during each setup. Examination of the GECDOE frequency 
plots supports this assumption. 

The results for frequency drift are shown in table U, and 
have been graphed in figure 1. They can be summarized as 
follows : 

a) The CM A -75 1 had a fairly uniform value for frequency 
drift, using either satellite, and easily met its spa- 
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TABLE 4 LONG-TERM OSCILLATOR FREQUENCY DRIFT 1 


SETUP 

NO. 

SATELLITE 

NO. 

CMA-751 

JMR-1A n 

JMR-1A n 

MX-1502 

1 

14 

0.45 

1.65 

3.15 

0.41 2 3 


19 

0.23 

1.78 

3.22 

0.11 2 

2 

14 

0.14 


2.50 2 

3 


19 

-0.32 


2.66 2 

3 

3 

14 

0.50 


2.57 2 

-2.88 


19 

0.33 


2,87 2 

-2.72 2 

4 

14 

0.39 


2.27 2 

0,47 


19 

0.10 


2.68 2 

1.11 

5 

14 

0.27 


2,06 

0,78 2 


19 

0.50 


2.47 

0.74 2 


SPECIFICATION; 


/ DAY 

±1.00 

±0,50 

±0.50 

? 

/100 S 

±0.01 

±0.05 

±0.05 

±0,08 


1 10“ 10 PARTS PER DAY. DETERMINED FROM FREQUENCY OFFSET OF FIRST 
AND LAST PASS OF SINGLE STATION. PRECISE EPHEMERI S SOLUTION. 

2 SOLUTION SHOWS FREQUENCY JUMP AFTER FIRST OR SECOND PASS. 

3 TOO FEW PASSES IN SOLUTION. WITH TWO FREQUENCY JUMPS 
(OSCILLATOR DISTURBED DUE TO MAINTENANCE) 
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Fig. 4 Oscillator frequency drift versus time. 
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cified IQ” 10 parts/day precision. Ihe frequency drift 
uas usually from one half to one tenth of that value, 
and even approached its 100 second specification, 

b) 'Ihe JMR-1A #2 also had a fairly uniform value for fre- 
quency drift, using either satellite. However, both 
it and the JKE-1A <M failed to meet their 0.5 x 10“i° 
parts/day specified precision, (Note that this speci- 
fication is actually fcr the JMB-1. It is assumed 
that the J ME- 1 A would have the same or a fetter speci- 
fication. ) 

c) Ihe MX1502 did not have a consistent value fcr its 
frequency drift, which shews oscillations during the 
second through fourth setups. Since the values for 
the first, fourth and last setups are at least simi- 
lar, one would suspect that the frequency drift chang- 
es are mostly due to the various times that the in- 
strument was opened (and its oscillator turned off) 
for repairs. No specification for the MX1502 drift 
per day is available for comparison purposes. 


IINAI_CCHMENTS 

The results presented here should be considered as the 
final ones of this comparison, although if time permits, 
some additional material will possibly be added to the re- 
port version of £Archinal, 1982] and the final version of 
£ Archinal and Mueller, 1982]. Work is also continuing on 
the documentation and further testing of the IEh version of 
the GECDCP Program System. 

As to the further use of the data Detained, the recommen- 
dation is made here that the data from both this comparison 
and the Ottawa comparison be finally processed together in 
multi-station solutions, to provide a comparison cf hew well 
the various possible instrument pairs can measure the long 
Columbus-Gttawa baselines involved. furtner, it is also 
suggested that a similar reduction be made (if the data can 
be obtained) using the "Quebec" data described in [Rcceau, 
1981], which was also obtained during the operational phase 
of this comparison. 

Other investigations are also possible, including extend- 
ing the results given above by making further comparisons of 
the chords, comparing the vertical and horizontal positions 
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of the stations thcough theic cooLdinates, and ccmpariug the 
computed coordinates with the available ccntccl coordinates. 
These items were not done in this study mainly because they 
ace considered to be of lesser importance than the other re- 
sults presented# and due to a general lack of time ter these 
lengthy investigations. Other work concerning program op- 
tions or comparisons o£ programs could also be dene with 
this data. 
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2.3 Earth Deformation Considerations for the Maintenance of a 
Conventional Terrestrial Reference System 


The role of deformation analysis in the maintenance of a new Conven- 
tional Terrestrial Reference Frame has been outlined in previous semiannual 
reports and in [Bock and Zhu, 1982]. Basically, a set of fundamental 
coordinates x 0 of a global network of stations adopted at an initial epoch 
define the reference frame. The initial size and shape of this network 
is defined by the corresponding baseline lengths, D 0 . By comparing uhe 
estimated baseline lengths at a later epoch to Do, the deformations of 
the network can be estimated. This information is then used to improve 
the global estimates of variations in polar motion and earth rotation, 
with respect to the conventional axes defined by X 0 . 

Mathematical Model and Preliminary Estimation Model 

The mathematical model for the deformation analysis is simply the 
chord length of baseline i-j 

D u = c < x rV 2 + (y r y i )2 + (z j- z i )2]1 

This model is linearized about X 0 to yield 
L = AX + V 

1L 

where the observation vector L for the k tn baseline is 


and the parameter vector X represents the deformations, i.e., the change in 
coordinates between the initial epoch and a later one. V denotes the noise 
vector. 

Since the design matrix A is rank deficient by 6* we are restricted 
to a Generalized Gauss-Markoff (GSM) model. (L, AX, oo 2 ?" 1 ) where 

E(L) = AX 
D(L) = o 0 2 P _1 
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If there Is no a priori Information for the deformations, the minimum 
bias P-least squares estimate for X is given by 

X = N + U = Apj L-, N * A T PA, U = A T PL 

using the notation of [Rao and Mitra, 1971], where P is the weight matrix 
of the observations. This estimate can be shown to be equivalent to that 
obtained from augmenting the normal matrix N by a set of constraints C 
such that [Blaha, 1971] 

AC T = 0 
CX = 0 
and 

X = (N + C T C) -1 - C T (CC T CC T ) _1 C 

This means that we constrain the origin and orientation defined by the 
coordinates at some later epoch t to be equivalent to that defined by X 0 . 

Extended Models for A Priori Deformation Information 

In the case of the availability of a priori information on the 
deformations of the network, e.g., as provided by absolute plate motion 
models, four possible estimators have been outlined and analyzed in 
[Bock, in preparation]. We briefly outline here the corresponding estimates 
and their respective properties. 

Consider an expanded GGM model (L, AX, Q v , Q-) where 

E(L) » AX 

D(L) = Q v = E£VV T > = a 0 2 P _1 
E(X X T ) = Qx 

* Zj + VJ (jijjr = E{ X } = X) 

where J is an independent estimate of the parameter vector. The resulting 
minimum M-norm P-least squares minimum variance estimate for X 

X x = Q X N(NQ X N) + U 
= M" 1 N(NM _1 N) + U 
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where M = Q x ~ x (positive definite). Xi has the property of minimum bias. 
Therefore, this estimate is termed the BLIMBE (Best Linear Minimum Bias 
Estimate). In this case, it can be shown that this estimate is equivalent 
to that obtained from augmenting the normal equations by CM such that 

AC T 3 0 
CMXi= 0 
and 

Xi 3 [(N + MC T CM)" X - C T (CMC T CMC 1 )" 1 C]‘ 1 U 

Therefore, we can say that the reference frame is maintained in a minimum 
M-norm P-least squares sense by a specified number of CTS stations. 

For positive semidefinite Q x , which would be the case for any plate 

model 

Xi 3 (N + M)" 1 N[N(N + M)' 1 N] + U 

with M = Q x + , In this case, the estimate is minimum M-seminorm P-least 
squares but is no longer minimum bias. 

For the BLIMBE we assume that the parameter vector X is deterministic 
and define a weighted norm in the parameter space on the basis of a priori 
information on X, Another possible biased estimator can be obtained by 
considering X as a random variable. Our estimation model is (L, AX’, Qy, Q x ) 
where 

E{ X > = J 

D[X] = E{(X-X)(X-X) T } 3 z x 

which gives 

Q x = E{XX T > = s x + X X T 

The vector X includes the deformations computed from, say, a plate motion 
model and z x is its covariance matrix. The distribution of L is given by 

E{L } 3 AX 

D[L] 3 A E x A T + o 0 2 P _1 

from which 

Q L 3 E{LL T } 3 A Q x A T + a 0 2 P‘ x 
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Q xl = E{XL T } - Q x A T 

and we assume 

Q xv = 0 | 

By the Gauss-Markoff theorem [Liebelt, 1967] * 

* 2 - Q XL q L ‘‘ L 

■ Q x A T (A Q x A T + a, 2 ?" 1 )' 1 L 
Which, for positive definite Q x , 

X 2 = (N + M)" 1 U : M = Q x “ x 

This estimate has been referred to as the Best (or Bayes) Linear Estimate 
or BLE for short [Rao, 1973, 1976]. While the BLIMBE has the minimum bias 
property, the BLg has minimum mean square error, i.e., it minimizes the 
sum of covariance and biased squares 

MSE(X) = E x + [X - E(X)][X - E(X)] T | 

in the class of biased estimators. Note that the BLE requires some jj 

knowledge of the deformations in order to compute Q x . Furthermore, while 
the BLIMBE reference system is maintained through the constraints CM Xl = 0, j 

the deformations estimated by the BLE are with respect to an underlying jj 

reference frame of *he deformation model from which Q x is computed. \ 

The previous two estimates are drawn from the class of biased esti- 
mators. If Q x can be constructed, that is, if there exists a priori defor- 
mation information, then the origin and orientation singularities are 
essentially eliminated. We then are led to investigate whether an unbiased 
estimate exists and we find the Bayesian estimate. Consider the estimation 
model (L, AX, Qy, X, Zj) where X is deterministic, X random and the set of I* 

observation equations 

' L 1 A! [V " 

= X + ; L y = X 

N M W 

such that 
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E{X> = X 

D[X] = E{(X - X)(X - X) T > = Zj 
E{L} = AX 

D[L] = A Zj A T + a 0 2 P _1 = Z L 

The least, squares solution for this model yields 

X 3 = Zj A T (A Zj A T + El)’ 1 L 

+ [I - Zj A T (A e^ A T + e l )’ x A] X 

for E^ positive semidefinite. For positive definite E-^» this reduces to 

X 3 = (N+M)" 1 U + [I - (N+M)" 1 N]X ; M = zf l 
= X + (N+M)’ 1 A T P(L - AX) 

= (N+M)’ 1 (U + MX) 


It is easiliy seen that given this estimation model, particularly 
E{X} = X, E{X 3 } = X, so that X is unbiased. This estimate has the minimum 
mean square error property which implies minimum variance since the bias 
is equal to zero. Note that in the BLE, the a priori information is incor- 
porated i.lto the moment matrix Q x , while for X 3 , "X is applied directly, and 
a residual deformation is estimated. Thus, we can consider the BLE (X.a) 
as a "weak." Bayesian estimate and X 3 a "strong" Bayesian estimate. 

Assume again that some a priori deformations are available, In this 
case, the model may indicate that CX=L X where f 0 which leads to an 
alternative approach to the constraints CMX x = 0 of BLIMBE. Consider the 
following set of observation equations 

L x = CX 

We assume the estimation model (L,AX | CX = CX, e x , Q v ) where 

Ei'CX} = CX 

D[CX] = C Zj C T 

E{L} = AX 

D[L] = A Zj A T + oo 2 ?' 1 
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For this model, the least squares estimate Is 


X H » [N + C T P X C]" 1 U + C T P x C 1 

where 

p x = (C Z-^ C T )” 1 


From [Chipman, 1964] 

Ap| * [N + C T P x C]" 1 A T P 

Cp j = [N + C T P x C]" 1 C T P x 
X 

so that 


X 4 n Apj L + Cp j X 

A 


i + u + cj T J 
p x* 


Therefore, X 4 can be viewed as a correction term to the minimum I-norm 

* 

P-least squares estimate X*, or a combination of the BLIMBE and Bayesian 
approaches. 


The properties of the four estimators are summarized in Table 1. 


Addition and Temporary Deletion of CTS Stations 

The reference frame is defined by a particular number of CTS stations. 
It is quite possible that from time to time one or more of the stations will 
not be able to participate in a particular deformation analysis observing 
session which should involve all stations. Furthermore, it must be anti- 
cipated that new stations will be added to the frame periodically. Both 
of these occurrences must be dealt with in order to maintain continuity 
and avoid ambiguity in the reference frame definition. For the addition 
of CTS stations we use the filtering and estimation capabilities of least 
squares collocation. The model becomes 

L = AX + BS + V 

Where X is deterministic and represents the coordinates of the new stations 
to be estimated. The vector S, the signal, is random and includes the 
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filtered deformations. The L and V vectors are as before. From 
[Moritz, 1980), 

X = [A T (BQ s B T + Q V PA] + A(BQ s B T + Q„P L 
S = q s B T (BQ s B T + Q V P (L - AX) 

where Q s is the same as the previous Q x * 

If a station cannot observe, we can use the prediction capabilities 
of least squares collocation to precict tne deformation via 

S = Q s B T (BQ t B T + Q v ) _1 L 

where 

~t ' 

s= U- 

t includes the deformation of the observing station, and u the predicted 
deformations of the missing stations. 

Conclusions 

In order to test the properties of the four estimators and their 
suitability in estimating deformations, a series of simulations were run 
as described in [Bock, in preparation]. A 20-station, 8-plate network was 
chosen for the simulations as depicted in Fig. 1 and Table 2. The AMI-2 
absolute plate motion model of [Minster and Jordan, 1978] was "adopted" 

(see Table 3) . 

The following conclusions were arrived at based on the simulations. 
Assuming that the absolute motion models available today are good to within 
their stated noise levels (this is reasonable considering that [Bender, 

1981] indicates that their predicted deformations differ at the centimeter 
level), it is found that it is advisable to adopt a deformation model than 
not at all. This was seen from comparing the deformation estimates obtained 
with a deformation model and those obtained when M = I (no model) is assumed. 
If a model is adopted, then the BLE appears to be the best candidate for 
deformation analysis. This conclusion follows from several considerations. 
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Table 2 20-Station, 8-Plate Simulation Network and OF POOR QUALITY 
AMI-2 Velocities 


No. 

Station 

Latitude 
D M 

Longitude 
D M 

Plate 

Velocity (cm/yr) 
X x Z 

1 

Johannesberg* 

-29 

10 

28 

02 

AFBC 

0.03 

0.99 

1.00 

2 

Cairo 

30 

3 

31 

15 

A PRC 

-0.49 

-0.31 

1.01 

3 

Xagos 

6 

27 

3 

23 

AFRC 

-0.09 

0.34 

0.62 

4 

Onsala* 

57 

0 

13 

0 

EU8A 

-0. 14 

-0..32 

0. 14 

5 

Jodrell Bank 

53 

0 

358 

0 

EUR A 

-0. 13 

-0.31 

0.09 

6 

Snanqhai 

30 

45 

121 

45 

EURA 

-0.09 

-0..20 

0.21 

7 

Fairbanks 

64 

50 

212 

10 

NOAM 

-1.38 

-0.15 

-0.59 

8 

Ft. Davis* 

30 

38 

256 

3 

NOAM 

-2.43 

-0.07 

-1.11 

9 

flestf ord 

42 

36 

288 

30 

NOAM 

-2.27 

-1.29 

-0.54 

10 

Maui* 

20 

30 

203 

45 

PCFC 

-1.80 

8.35 

4.59 

11 

Tahiti 

-17 

30 

210 

30 

PC FC 

-6. 11 

7,5 9 

4-48 

12 

Marshall Isles 

7 

0 

167 

0 

PCFC 

2.7 3 

9.24 

4.76 

13 

Sao Paulo* 

-23 

33 

313 

21 

SO AM 

-2. 2b 

-1.94 

-0.33 

14 

Buenos Aires 

-34 

37 

301 

36 

SOAM 

-2.44 

-1.30 

-0.25 

15 

Caracas 

10 

35 

293 

40 

SOAM 

-2.75 

-1.26 

-0.26 

16 

Orroral* 

-35 

18 

149 

8 

INDI 

-3.. 63 

1.70 

5.63 

17 

Yaragadee 

-29 

3 

115 

21 

INDI 

-4.20 

1.99 

6.47 

18 

Bombay 

13 

56 

72 

51 

IND.I 

-0„ 95 

-1.25 

4.30 

19 

Easter Isle* 

-27 

5 

250 

39 

NAZC 

6., 04 

-1.56 

-1.04 

20 

| 

Arabia* 

24 

39 

46 

46 

ARAB 

-1.42 

-0.36 

2.70 


* 8-Station 8- 

Plate 

Network 









Table 3 AMI-2 Absolute Motion Plate Model (Adapted from 
[Minster and Jordan, 1981], Table 7) 




Absolute Rotation Vector 


Plate 

Deg 

<N) 

Deg 

(E) 

Deg/M- 

Y. * 

1. African 

18.76 

33.93 

338.24 

42.20 

0. 139 

0.055 

2. Eurasian 

0.70 

124.35 

336,81 

146.67 

0,038 

0.057 

3. North American 

-50.31 

16.21 

319. 33 

39.62 

0.247 

0.030 

4. Pacific 

-61.66 

5. 11 

97. 19 

7.71 

0.967 

0.085 

5- South American 

-82.28 

19.27 

75.67 

85.88 

0. 285 

0.084 

6. Indian- Australian 

19.23 

6.96 

35.64 

6.57 

0.716 

0.076 

7- Nasca 

47.99 

9.36 

266.19 

8. 14 

0.585 

0.097 

8. Arabian 

27.29 

12.40 

356.06 

18.22 

0. 388 

0,067 

9. Antarctic** 

21-85 

91.81 

75.55 

63.20 

0. 054 

0.091 

10. Caxribean** 

-42.80 

39.20 

66.75 

40,98 

0. 129 

0. V34 

11. Cocos** 

21.69 

3.08 

244.29 

2-81 

1.422 

0. 1 19 

* Million Years 





. 


** Not used in the 

Simulations 








First, the BLE provides the best estimates (in the sense of minimizing 
the root mean square error between true and estimated deformations) at 
the same level as the Bayesian estimate, in the case when the deformation 
model is correct (and then the deformation is just being filtered from the 
baseline noise). Second, and most important, it is markedly less sensitive 
to errors in the adopted deformation model. This is particularly 
apparent in the case that in reality there is no deformation but we assume 
some deformation model. These results are due to the minimum mean square 
error and minimum norm properties of the BLE and its "weak" Bayesian inter- 
pretation. 

Finally, we should stress that the reference system is dependent on 
the choice of estimation models including the choice of M (as well as P, 
but to a lesser extent). This leads to the need for investigations 
concerning how sensitive the reference system is to changes in M and P. 

For example, what measures should be taken as M and P improve. with time. 

The algorithms presented here are general enough to incorporate geo- 
physical as well as geodetic evidence of deformations. In [Bock, >n prep- 
aration] only models for deformations of interplate type have been con- 
sidered, to be monitored by periodic re-observations of the baseline lengths. 
Other aspects to be considered include intraplate and local motions (the 
site stability problem). Local effects can possibly be modeled on the 
basis of on-site observations such as by tidal gravimeters and local 
geodetic nets. It is necessary to investigate how to incorporate these 
and other types of observations (and their corresponding reference frames) 
into CTS operations. 

This investigation is now being completed, and the final report is 
in preparation by Y. Bock, to appear in the report series of the Depart- 
ment of Geodetic Science and Surveying, The Ohio State University. 
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2.4 Development of Models for Studying Ice Sheet and Crustal 
Deformations 


The observed locations of survey markers change with time. When 
random and systematic errors are accounted for, what remains is actual 
movement. The movements of a network of stations can be described as the 
translation and rotation of the stations as a group and the deformation 
occurring within the network. Thus when a network of stations is resurveyed, 
it should be possible to obtain the geophysical parameters of velocity, 
rotation rate and strain rate [Dermanis, 1981; Livieratos, 1980; Reilly, 
1979]. If the same network is resurveyed more than once, either the 
derivatives of these quantities or averaged values may be calculated. 

As most stations are on the surface of the earth, it is natural to 
assume that all movements and deformations are two-dimensional. This may 
oe adequate in many cases. However, vertical movement and deformation may 
occur because of irregularities in the surface, faulting, or from being 
buried under new material. Also, for networks covering relatively large 
areas, the surface of the earth cannot be well approximated by a plane. 

In this case, it may be better to determine the movements in an arbitrary 
(earth-centered) coordinate system and then transform these results to a 
latitude, longitude and elevation coordinate system. 

A model is being developed to determine these geophysical parameters 
from the coordinates of a network that has been resurveyed at least once. 
Several methods have been proposed for obtaining sufficiently accurate 
coordinates [Brunner et al., 1981, Niemeier, 1979]. One technique that 
has been proposed for studying tectonic deformation is to use positions 
determined by Doppler satellite receivers [Malyevac and Anderle, 1979]. 

The precision of the receivers used individually (point positioning) is 
meters to tens of meters. But by using translocation between two or more 
receivers, the relative positions can be determined to within decimeters 
[Brown, 1976]. However, the movements and deformations of the crust are 
slow even in tectonically active areas [Savage, 1978; Minster and Jordan, 
1978]; thus the time span between resurveying must be of the order of 
decades. Because the time period between reobservations is so long, it 
may be difficult to guarantee that the coordinate systems are identical. 
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For example, the coordinate system defining the broadcast ephemeris of 
the Navy Navigational Satellite System slowly varies with time, This 
problem could be overcome by using relative rather than absolute coordinates. 
Thus the velocities and rotation rates would be relative to some "fixed" 
stations. However, the deformation within the network can still be obtained 
by calculating the strains from the changes in the chord lengths between 
the stations. The only assumption needed for this is that the scale of the 
coordinate system has not changed, Because the strains obtained this way 
are theoretically identical to the strains obtained from coordinate differ- 
ences, any differences can be attributed to rotations and/or translations 
of the coordinate system. 

For the purposes of testing the model, the data set being used is 
from survey stations placed on the Greenland ice sheet. Seven Magnavox 
1502 satellite receivers were used during the summers of 1980 and 1981 to 
obtain the movement of 22 stations on the ice sheet of Greenland. Using 
the data reduction program GEODOP [Kouba and Boal , 1976], the coordinates 
of the stations have been obtained relative to the positions of two 
stationary stations (which were located on the west coast of Greenland). 

The formal accuracy of the coordinates is under 20 cm. These stations 
are moving at velocities of up to 45 m per year, and the magnitude of the 
maximum strain rates are over 100 ppm. 
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